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ABSTRACT 

We present new measurements of the evolution of the X-ray luminosity functions (XLFs) 
of unabsorbed and absorbed Active Galactic Nuclei (AGNs) out to z ~ 5. We construct 
samples containing 2957 sources detected at hard (2-7 keV) X-ray energies and 4351 sources 
detected at soft (0.5-2 keV) energies from a compilation of Chandra surveys supplemented 
by wide-area surveys from ASCA and ROSAT. We consider the hard and soft X-ray samples 
separately and find that the XLF based on either (initially neglecting absorption effects) is best 
described by a new flexible model parametrization where the break luminosity, normalization 
and faint-end slope all evolve with redshift. We then incorporate absorption effects, separately 
modelling the evolution of the XLFs of unabsorbed (20 < log Wh < 22) and absorbed (22 < 
log Nn < 24) AGNs, seeking a model that can reconcile both the hard- and soft-band samples. 
We find that the absorbed AGN XLF has a lower break luminosity, a higher normalization, 
and a steeper faint-end slope than the unabsorbed AGN XLF out to z ^ 2. Hence, absorbed 
AGNs dominate at low luminosities, with the absorbed fraction falling rapidly as luminosity 
increases. Both XLFs undergo strong luminosity evolution which shifts the transition in the 
absorbed fraction to higher luminosities at higher redshifts. The evolution in the shape of the 
total XLF is primarily driven by the changing mix of unabsorbed and absorbed populations. 

Key words: galaxies: active - galaxies: evolution - galaxies: luminosity function, mass func¬ 
tion - X-rays: galaxies 


1 INTRODUCTION 

The luminosity function of Active Galactic Nuclei (AGNs) repre¬ 
sents one of the crucial observational constraints on the growth of 
supermassive black hole (SMBHs) over the history of the Universe. 
The shape of the luminosity function reflects a combination of the 
underlying distribution of the SMBH masses a nd the distribution 
of their accretion rates or Eddington rat ios (e.g. lAird et al.ll^l3al : 
IShankar et al.ll2013l : ISchulze et al.ll2015h . Thus, accurate measure¬ 
ments of the shape and evolution of the luminosity function provide 
crucial insights into the physical processes that drive SMBH growth 
over cosmic time. 

Many AGNs are surrounded by gas and dust that can obscure 
their emission at certain wavelengths. Thus, it is vital to under¬ 
stand AGN obscuration in order to obtain accurate measurements of 
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the luminosity function. Quantifying AGN obscuration also reveals 
whether SMBHs undergo significant periods of obscured growth, 
when this takes place within the lifetimes of AGNs, and how it re¬ 
lates to the triggering and fuelling processes. 

Obtaining accurate measurements of the luminosity function 
and revealing the extent of obscuration requires large, unbiased 
samples of AGNs selected over the widest possible range of red- 
shifts and luminosities. Optical surveys, combined with follow- 
up spectroscopy , can efficiently cover wide areas (e.g. SDSS: 
lYork et aTIbOOol) but are biased towards the most luminous, unob¬ 
scured sources. Alternatively, AGNs can be identified in the mid- 
infrared (mid-IR), which probes the reprocessed emission from the 
dusty, obscuring material. Mid-IR selection should not be biased 
against obscured sources, but the contribution of the host galaxy is 
often significant at these wavelengths, which limits mid-IR selec¬ 
tion to lumino us sources where the galaxy light is over whelmed by 
the AGN ie.g. lDonlev et al.l2008l : lMendez et al.l2013h . 
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X-ray surveys can efficiently identify AGNs over a wide 
luminosity range, including low-luminosity sources where the 
host galaxy domi nates at optical or infrared wavelengths (e.g. 
iBarger et al.ll20()3l) . Nevertheless, soft X-ray emission (at energies 
< 2 keV) will be absorbed by the same gas and dust that obscures 
the AGN at optical and UV wavelengths. Thus, soft X-ray sam¬ 
ples are generally dominated by unobscured AGNs. Absorption bi¬ 
ases are reduced at hard X-ray energies 2 — 10 keV), except in 
the most heavily-obscured, Compton-thick AGNs (equivalent line- 
of-sight hydrogen column densities A^h ^ 10^"^ cm“^). However, 
even Compton-thick sources may still be identified at soft or hard 
X-ray energies due to scattered emission, including the Compton- 
scattered emission (“reflection”) from the obscuring material itself. 

A large number of deep and wide X-ray surveys have been car¬ 
ried out, taking advantage of t he efficiency and power of X-ray se¬ 
lection (see a recent review bv ISrandt & Alexandeill201^ . A num¬ 
ber of studies have measured the X-ray luminosity fun ction (XLF) 
of AGNs out to high redshifts using these s amples (e.g. lUeda et al.l 
l2003l : [Barger et al.T2005l : lMivaii et al.l2015h . These studies find that 
AGNs are a strongly evolving population, with a sharp decrease in 
their number density between 2 ; ~ 1 — 2 and today. Bright X-ray- 
selected AGNs are found to peak in number density at 2 « 2, sim¬ 
ilar to optically-selected QSOs. Fainter AGNs peak later in the his¬ 
tory of the Univ erse (2 ~ 1) but with a much milder decline to the 
present day (e.g. lHasinger et alj2005h . These patterns have led sev¬ 
eral authors to propose a luminosity-dependent density evolution 
(LDDE) pa rametrization to describe the evolut ion of the XLF of 
AGNs (e.g. lMivaii et al.|[200d : IUeda et alj|2003h . In this model the 
XLF is modified by differing degrees of density evolution that vary 
with luminosity and redshift. This results in an XLF that changes 
shape o ver cosmic time. 

In lAirdetalJ 1201(1 hereafter A10) we challenged this evo¬ 
lutionary model with a detailed study of the hard-band XLF that 
carefully accounted for numerous uncertainties and biases that 
were generally not included in prior measurements. These included 
flux measurement errors, Eddington bias, incompleteness of opti¬ 
cal identifications, and the uncertainty in photometric redshift esti¬ 
mates. At high redshifts (2 ~ 2 — 3) we adopte d a rest-frame UV 
colour pre-selection technique dAird et alj|2008h . By performing a 
robust model comparison based on Bayesian statistical techniques, 
we found that the evolution of the XLF could be described by a sim¬ 
pler model in which the XLF retains the same shape at all redshifts 
but evolves in both lu minosity and density (see also lAssef et al.l 
|201l|;lRosse7^l2013h . 

While A10 presented a number of important advances, ab¬ 
sorption effects were not explored. Other studies have attempted to 
measure the distribution of absorption column d ensities and presen t 
absorption-corrected measurements of the XLF. IUeda et alj l l2003ll 
found that the fraction of absorbed AGNs (those with Ah > 10^^ 
cm“^) was strongly dependent on luminosity, decreasing at higher 
luminosities. Later studies found that the fraction of absorbed 
AGNs depends on both luminosity and redshift, dropping at high 
luminositie s but increasing (at a given luminosity ) to higher red¬ 
shifts (e.g. IlT Franca et all l200^ : lHasing^ 120081) . The extent of 
any redshift evolution has be en a matter of debate (see lAkvlas et al.l 
120061 : lOwellv & Pa^bOOd) . potentially due to difficulties regard¬ 
ing the selectio n functions for ab sorbed and unabsorbed sources. 
Recent work bv lUeda et alj ^2014) re-examined the evolution of the 
XLF and the distribution of Ah (the “Ah -function”) using a large 
compilation of both soft and hard X-ray surveys and found that 
both a luminosity and redshift dependence of the absorbed frac¬ 
tion were required. They also found that an LDDE parametriza¬ 


tion was needed to describe the evolution of the XLF (with some 
further modifications to describe the evolution at 2 > 3, see also 
ICivano et alj20li] : lHiroi et al.ll201 j) . 

Recently, the combination of extremely deep X-ray survey 
data and new approaches to X-ray spectral analysis have en¬ 
abled improved measurements of Ah at 2 ~ 0.5 — 2 and 
have been u sed to identify sizab l e samples of Compton-thick 


AGNs (e.g. IComastri et al.l 


ty sizab l e samples ot Lompton-micK 
, _ lI 1 201 ll : iGeorgantopoulos et al. 201^ 

iBrightman et alj |2014|) . Ruilding on this work. IPuchner et al.l 
( I2OI5I) used a flexible, non-parametric method to estimate the 
space densities of AGNs as a function of redshift, luminosity and 
Ah, effectively measuring the XLF for different column densi¬ 
ties. This work also recovered a luminosity and redshift depen¬ 
dence in the evolution of the fraction of absorbed AGNs (although 
the Compton-thick fraction was consistent with a constant value 
of 35 per cent). However, a detailed comparison of parametric 
models for the evolution of the XLF of AGNs was not undertaken. 

In this paper we address some remaining issues in studies of 
the evolution of the XLF of AGNs: the shape of the XLF and how 
it evolves with redshift, the extent of any luminosity and redshift 
dependence of the absorbed fraction, and the connection between 
the absorption properties and the evolution of the AGN population. 
We combine samples selected at both hard and soft X-ray energies 
and determine the underlying XLF and distribution of Ah that ad¬ 
equately describ es the observed fl uxes in both samples (similar to 
the ap proach of lUeda et alJ l2014L cf. the X-ray spectral analysis 
used in iBuchner et alj2015ll . 

In Section 1^ we describe our datasets that we use to define 
large samples of X-ray sources selected in the hard (2-7 keV) and 
soft (0.5-2 keV) energy bands. We also compile deep optical, near- 
infrared and mid-infrared imaging across our fields that we use to 
robustly identify counterparts to our X-ray sources and calculate 
photometric redshifts. In Section[3 we describe our Bayesian sta¬ 
tistical technique that allows us to incorporate a range of X-ray 
spectral shapes and account for the effects of absorption. We also 
introduce an approach to account for the contribution from nor¬ 
mal, X-ray detected galaxies on our measurements. We then present 
measurements of the XLF based on our hard and soft samples in¬ 
dividually (Section |4}, introducing a new flexible parametrization 
of the XLF. We show that significant discrepancies between the 
measurements at all redshifts warrant the further consideration of 
absorption effects. In Section [5] we separately model the XLF of 
unabsorbed and absorbed AGNs (including a contribution from 
Compton-thick sources) and show how the combination of these 
populations can simulataneously account for both our hard- and 
soft-band samples. Our results place constraints on the total XLF 
of AGNs and the absorbed fraction as a function of luminosity and 
redshift. In Section|^we compare our results to prior work and dis¬ 
cuss the wider implications of our findings. Section|7]summarizes 
our paper and overall conclusions. 

Given the length of this paper, a casual reader may wish to skip 
to Section l4!4l Section|5]and the discussion in Section|^(and focus 
on Figures [T] [8] an d[9). We adopt a flat cosmology with Ha = 0.7 
and /i = 0.7 throughout this paper. 


2 DATA 

To constrain the evolution of the XLF we require large samples of 
X-ray selected AGNs. By selecting samples in both the hard (> 
2 keV) and soft (0.5-2 keV) observed energy bands, we can also 
constrain the distribution of Ah and correct for these effects on the 
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XLF. In this paper we combine a large number of Chandra X-ray 
surveys along with larger area surveys from A5CA and ROSAT. We 
give further details of our datasets below. Table [T] summarizes the 
different surveys and provides the number of hard and soft X-ray 
selected sources from each. 


2.1 Chandra X-ray data 


In this paper we use Chandra X-ray observations from five dis¬ 
tinct parts of the sky: the Chandra Deep Field-South (CDFS), Chan¬ 
dra Deep Field-North (CDFN), Extended Groth Strip (EGS), COS¬ 
MOS, and Bootes fields. 

In the CDFS field we identify two different “surveys”: 1) the 
series of observations that have targeted the central ~ 0.07 deg^ of 
the field and reach a total combine d exposure time of ~4Ms (the 
CDFS-4MS survev. IXue et al.l201 ih : and 2) the series of four 250ks 
observations that s urround the central area (the Extended-CDFS, 
or E-CDES survey, iLehmer et al.ll2005ll . The X-ray data reduction 
and source detection procedures were carried out independently for 
each of these surveys. To combine the two surveys, we define a cen¬ 
tral region where the CDES-4Ms survey takes precedence, roughly 
corresponding to the area within ~ 9' of the center of the field. 
Outside this central region we adopt the sources detected in the E- 
CDFS survey only. This procedure ensures we have a well-defined 
sample where we can accurately determine the sensitivity. 

In the EGS field we also identify two distinct surveys: 1) the 
series of eight pointin gs that reach a no minal depth of ~200ks 
and were presented in iLaird et ^ ( l2009l) . which we refer to as 
the AEGIS-X W(ide ) survey; and 2) the AEGIS-XD(eep) survey 
dNandra et al.ll201^ which took three of the original 200ks to a 
depth of ~ SOOks. We adopt the deeper AEGIS-XD observations 
when available. 

In each of the remaining fields we adopt data from a sin¬ 
gle C handra survey: the 2Ms survey in CDF N ([Alexander et al.l 
2003h: the ^16Qks C-COSMOS observations ( lElvis et al.ll2009l: 


Puccetti et al.ll200^ : and the 5ks XBootes survey jMurrav et d] 


2005h . 


The X-ray data from all of our surveys were reduc ed with our 
own p ipeli ne procedure, which is described in detail bv iLaird et al.l 
l l2009h and iNandra et al.l ( l2015h . We pe rformed point sour ce detec¬ 
tion using the procedure described by iLaird et ^ ( l2009h and ap¬ 
plied a false Poisson probability threshold of < 4 x 10“® to gen¬ 
erate catalogues of detected sources in the soft (0.5-2 keV), hard 
(2-7 keV), full (0.5-7 keV) or ultrahard (4-7 keV) observed en¬ 
ergy bands. We combined the source lists in each band to create a 
merged catalogue, which is used in the counterpart identification 
procedure described in Section l2^ below. 

For the 5ks XBootes survey we applied a stricter false prob¬ 
ability cut (< 10“®) in addition to a requirement of ^ 5 total de¬ 
tected counts in each band. This cut reduces the sample size but ap¬ 
plies an effective X-ray flux limit that helps raise the completeness 
of the spectroscopic follow-up in this field. In this field, we also 
restrict our analysis to the ~ 7.1 deg^ of the Bootes field that cor- 
respon ds to the 15 standard s ub-fields of the AGES spectroscopic 
survey dKochanek et al.l2012l) . This cut also ensures we have a high 
spectroscopic completeness (see Section l23] below). 

We determined X-ray sensitivity maps and are a curves for 
each band as described in iGeorgakakis et al.l ( l2008l) . accounting 
for the stricter false probability cut and minimum counts require¬ 
ment for the XBootes survey. We convert the sensitivity maps to 
area curves as a function of flux by assuming a fixed X-ray spectral 
slope of r = 1.4. We note that the flux calculated with this fixed 


conversion factor scales directly with the count rate; in Section |3 
below we describe our procedure to convert between the count rate 
and intrinsic quantities (such as luminosity), which allows for a 
more complex X-ray spectrum and accounts for uncertainties in this 
conversion factor. The sensitivity map calculation is limited to the 
footprint of the multiwavelength photometry for each field. Eigure 
[T]shows the corresponding area curves for each of our X-ray fields. 

We note that previously published X-ray source catalogues, 
often including mutiwavelength counterpart information, are avail¬ 
able for the all of our Chandra fields ( e.g. [Alexander et al.ll2003l ; 
iGoulding et al.ll201^ [Brand et al.ll2006l) . Adopting our own X-ray 
reductions and source detection procedures ensures we can accu¬ 
rately determine the sensitivity in a consistent manner, which is 
essential for our Bayesian analysis of the XL E. Our catalogues 
con tain ^ 10 — 25 per c ent fewer sources than IXue et al.l ( 1201 ill 
and IPuccetti et al] l l2009ll in the CDFS-4Ms and C-COSMOS ar¬ 
eas, mainly due to our stricter (i.e. lower) false probability thresh¬ 
old. Thus, our catalogues are more conservative. 

2.2 Multiwavelength counterparts and photometry 

In each of our Chandra fields, we identify multiwavelength coun- 
terpa rts to our X-ray sour c es using the l ikelih o od ratio (LR) metho d 
(e.g. ICiliegi et aD I 2 OO 3 I: lBrusaetal.1 l2007l; ICivano et al.l |2012|) . 
matching to multiple optical, near-IR and mid-IR bands to ensure 
a high completeness and reliability. We also compile multiwave¬ 
length photometry from a larger number of bands, which we use to 
calculate photometric redshifts (see Section|2j6]below). 


2.2.1 CDFS, CDFN and EGS 

In three of our fields—the CDFS, CDFN and EGS—we identify 
counterparts and extract the multiwavelength photometry using a 
custom version of the Rainbow Co s mological Surveys Databasj^ 
dPerez-Gonzalez et al.l 120051 l2008l ; iBarro et ^ 1201 lal f9). which 
provides a compilation of the various photometric datasets. Ap¬ 
pendix lists all the different photometric imaging datasets for 
each field that are used in this paper. We note that full coverage of 
the entire field is not always available in each photometric band. 
All the images are registered to a common astrometric reference 
frame and photometry is performed in consistent apertures to pro¬ 
duce spectral energy distributions that span from the UV to mid- 
IR. The X-ray sour ce matching procedure is described in detail in 
iNandra et alj ( 12015h . but we briefly summarise the method here. 

Eirst, we reran the Rainbow photometric code, extracting all 
potential counterparts within 3.5” of the X-ray positions in any of 
the bands covered by Rainbow (using initial SExtractor catalogues 
in each of the bands). The counterparts were then cross-matched 
using a 2” search radius to create a single multiband catalogue. 
We obtained consistent photometry by applying a single aperture 
across all optical and near-IR bands. We als o extracted IRAC pho¬ 
tometr y, ap plying the procedure described in IPerez-Gonzalez et al.l 
(I2OO8I) and iBarro et alJ ( 1201 lal) to deblend the IRAC photometry 
when a single IRAC source is associated with multiple optical/near- 
IR counterparts. 

Next, we calculated the LR for all candidate counterparts de¬ 
tected in the IRAC 3.6pm band and determined an LR threshold 
that maximises the sum of the completeness and reliability (see 
ILuo et aLll201(]h . A candidate counterpart that exceeded this LR 


^ https://rainbowx.fis.ucm.es 
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Table 1. Details of the X-ray surveys used in this paper. 


Field 

(1) 

Survey 

(2) 

RA 

(J2000) 

(3) 

Dec 

(J2000) 

(4) 

X-ray 

exposure 

(5) 

Survey 

area 

(deg2) 

(6) 

Soft band 

Ax Actrprt 

(V) (8) 

^spec —z 

(9) 

Hard band 

-^X -^ctrprt 

(10) (11) 

^spec —Z 

(12) 

CDFS 

CDFS-4MS 

03:32:27.2 

-27:47:55 

4Ms 

0.075 

413 

397 (96.1%) 

240 (58%) 

283 

273 (97%) 

155 (55%) 

CDFS 

ECDFS 

03:32:27.2 

-27:47:55 

250ks 

0.181 

334 

328 (98.2%) 

115 (34%) 

273 

268 (98%) 

101 (37%) 

CDFN 

CDFN 

12:36:49.3 

-h62:13:19 

2Ms 

0.112 

384 

363 (94.5%) 

242 (63%) 

286 

273 (96%) 

176 (62%) 

EGS 

AEGIS-XD 

14:19:20.8 

-h52:50:03 

800ks 

0.260 

698 

673 (96.4%) 

295 (42%) 

552 

539 (98%) 

233 (42%) 

EGS 

AEGIS-XW 

14:17:15.0 

-h52:25:31 

200ks 

0.204 

334 

332 (99.4%) 

130 (39%) 

274 

274 (100%) 

110 (40%) 

COSMOS 

C-COSMOS 

10:00:20.3 

-h02: 11:20 

160ks 

0.984 

1213 

1195 (98.5%) 

694 (57%) 

889 

877 (99%) 

530 (60%) 

Bootes 

XBootes 

14:31:28.3 

-h34:28:07 

5ks 

7.124 

754 

744 (98.7%) 

566 (75%) 

257 

255 (99%) 

196 (76%) 


ALSS 

13:14:00 

-h3 1:30:00 


5.800 




34 

34 (100%) 

33 (97%) 


AMSS 




81.77 




109 

109 (100%) 

107 (98%) 


ROSAT 




20391 

221 

221 (100%) 

221 (100%) 









Total 

4351 

4253 

2503 

2957 

2902 

1641 


Column (1) name of the field; (2) name of X-ray survey programme within this field; (3,4) approximate center of the suiwey; (5) nominal X-ray exposure 
time; (6) total area covered by both the X-ray data and the optical or infrared imaging, excluding areas close to bright stars; (7) number of X-ray sources 
detected in the soft (0.5-2 keV) energy band; (8) number (and fraction) of the soft X-ray sources that are associated with a robust multiwavelength 
counterpart; (9) number (and fraction) of the soft X-ray sources that have a spectroscopic redshift; (10) number of X-ray sources detected in the hard (2-7 
keV) energy band; (11) number (and fraction) of the hard X-ray sources that are associated with a robust multiwavelength counterpart; (12) number (and 

fraction) of the hard X-ray sources that have a spectroscopic redshift. 


threshold was deemed a secure counterpart (taking the counterpart 
with the highest LR value in cases of > 1 secure candidate). We 
then repeated the entire LR matching process for the bands indi¬ 
cated in the table in Appendix retaining any additional secure 
counterparts identified in these bands. This procedure enables us 
to identify secure counterparts for a high fraction (> 90%) of the 
X-ray sources in the CDFS, CDFN and EGS fields (see Table [Til. 
The vast majority (~ 92%) of the secure counterparts were iden¬ 
tified in the (deblended) IRAC 3.6p.m catalogue. Matching to the 
additional bands allows us to identify counterparts when the IRAC 
candidate is faint, blended, or non-existent. No additional cross¬ 
matching is required as the full multiband photometry is provided 
through matched apertures for all sources in the Rainbow database. 


2.2.2 COSMOS 


In the COSMOS field, which is not currently included in the Rain¬ 
bow surveys database, we matched directly between our X-ray 
source lists and two multiwavelength catalogues: 1) the COSMOS 
Intermediate and Broad Band Photometry Catalogue 200fl which 
is based on detectio n in the deep Suba ru imaging of the en¬ 
tire COSMOS field dCapak et al.ll2007l); and 2) the S-COSMOS 
IRAC 3.6pm based catalogue jSanders et al.l l2007h . Unlike the 
Rainbow catalogues, the S-COSMOS IRAC 3.6pm catalogue has 
not been deblended. Thus, we first identified secure matches (using 
the LR method) from the higher-resolution Subaru A catalogues. 
We found secure counterparts for 1348 of the 1621 X-ray sources 
(83 per cent) in our overall C-COSMOS catalogue. Next, we ap¬ 
plied the LR method to match between the X-ray and S-COSMOS 
catalogues. We identified secure counterparts for an additional 178 
sources, taking our overall completeness to 94 per cent. 

Finally, we cross-matched between our secure positions and 
the original catalogues, again using the LR method to ensure only 


secure associations are considered. For 1499 of the 1526 secure 
counterparts, we end up with both a Subaru and S-COSMOS 
counterpart. For these sources we obtain photometry in up to 18 
broad-band filters spanning from the UV to near-IR, as well as 15 
intermediate- or narrow-band optical filters, from the Subaru i+ 
catalogue (see lllbert et alj200 ^: lMcCracken et al .l2010h . The pho¬ 
tometry was extracted in 3” diameter apertures from PSF-matched 
images. We also adopt IRAC photometry in the 3.6, 4.5, 5.8 and 
8.0pm bands from the S-COSMOS catalo gue, applying the aper¬ 
ture corrections given in lllbert et al.l ( l2009l) . 

In 27 cases we have an A counterpart but do not find an S- 
COSMOS counterpart, either due to the limited depth of the IRAC 
imaging or because the i'^ source is blended at the IRAC resolu¬ 
tion. For these sources we simply ignore the IRAC bands in our 
photometric redshift estimates (see Section l231 below). 

In 23 cases, we identify an S-COSMOS counterpart, but no 
Subaru A source. For these sources we extracted photometry in 3” 
diameter apertures at the S-COSMOS position in the PSF-matched 
images for all the UV to near-IR bands, generally obtaining only 
upper limits for these optically faint sources. 


2.2.3 Bootes 

For the Bootes field, we compiled mult iwavelength catalogue s from 
the NOAO Deep-Wide Field Survey djannua & Devl[l9^ DR3, 
SDSS DR9 dAhn et alfcoH), (3ALFX CRtB FT.AMINCOS Ex- 
tragalactic Su rvey l Elstonet^i] 20061) . and the Spitzer Deep-Wide 
Eield Survey dAshbv et alJl2009 ). We use the LR method to match 
our X-ray catalogues to the appropriate selection band for each of 
these surveys, assigning secure matches from the surveys in the 
order of priority indicated in the table in Appendix The vast 
majority of our secure matches are identified in the NDWES / 


^ http://irsa.ipac.caltech.edu/data/COSMOS/datasets.html 


® http://galex.stsci.edu/GR6/ 
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Figure 1. X-ray area curves (sensitive area versus X-ray flux) for each of 
our surveys in the soft and hard bands. The black line indicates the total 
area curve for our study. We limit our sensitivity analysis to the area falling 
within the footprint of our multiwavelength photometry, excluding areas 
close to bright stars. The ROSAT area curve includes the ROSAT bright sur¬ 
vey and SA-N survey (see Section l23t . The EGS area curve includes both 
the AEGIS-XD SOOks data and the additional area at 200ks from AEGIS- 
XW. The CDES area curve combines the CDFS 4Ms data with the flanking 
250ks data from the E-CDES. These area curves assume a single, fixed con¬ 
version factor between the X-ray count rate and flux, conesponding to an 
X-ray spectrum with a photon index F = 1.4 and Galactic absorption only; 
the effect of different X-ray spectral shapes—and the resulting uncertainties 
in the conversion factors—are accounted for in our Bayesian methodology 
described in Sectionj^ 


band. We identify secure counterparts for 95.8 per cent of our X- 
ray sources. Finally, we cross-matched between the original cata¬ 
logues and our secure counterpart positions, again applying the LR 
method, to identify common sources. The combined surveys pro¬ 
vide photometry in up to 17 different bands (see Appendix lAt. 

2.3 Large-area surveys 


bright end of the XLF requires samples of higher luminosity X-ray 
sources identified from larger-area surveys. We thus supplement 
our sample with sources from large-area surveys carried out with 
ASCA (in the hard band) and ROSAT (in the soft band). 

For our large-area hard-band sample, we incl ude the 34 
sourc es from the ASCA Large Sky Survey (ALSS: lUeda et ^ 
Il999t) . which covers a contiguous area of 5.8 deg^ near the 
north Galactic pole. We adopt the optical identifications from 
lAkivama et alj ( l2000h : 2 sources are optically identified as galaxy 
clusters, 1 is a star, 1 source remains unidentified and the remain¬ 
ing 30 are associated with AGNs, all of which have spectroscopic 
redshifts. We also s elect sources from the ASCA Medium Sensitiv¬ 
ity Survey (AMSS. lUeda et al]|200lh . which combines data from 
a large number of ASCA observations at high Galactic latitudes 
over an area ~ 82 deg^. We include sources from the AMSSn sub¬ 
sample, selected in th e hard (2-10 keV) ban d, with optical identifi¬ 
cations presented by lAkivama et al.l l l2003ll . The sample includes 
87 X-ray sources and has 100 per cent spectroscopic complete¬ 
ness. We include additional sources with from the AMSSs sub¬ 
sample (Ueda & Akiyama, private communication), which includes 
20 AGN; 2 sources in this sample remain unidentified. We adop t 
are a curves from the A LSS and AMSS from lAkivama et all i2003h 
and lUeda et alJ l l2003h respectively. 

For our soft ban d sample we comb i ne samples from the 
ROSAT bright survey jFischer et alj 1 19981 : ISchwone et al. 2000h 
and t he Selected-Area-North survey (SA-N: lAppenzeller et al.l 
Il998ll . removing duplicate sources. We include sources with sig¬ 
nificant detections in the 0.5-2 keV band and adopt the unabsorbed 
flux estimates (corrected for Galactic absorption). We cut our sam¬ 
ple at flux limits of /0.5- 2 kev > 3.6 X 10 erg s ^cm ^ for 


the SA-N sample. These high flux limits ensure our sources all lie 
well above the sensitivity limits of the surveys, allowing us to adopt 
simple sensitivity curves that correspond to the entire area of each 
survey and cut off sharply at each of the flux limits (see Figure [TJ. 
All the sources above these flux limits have spectroscopic classi¬ 
fications and we identify a total of 221 AGNs with spectroscopic 
redshifts (excluding BLLac type objects). 

We note that, in contrast to our Chandra fields, our area curves 
for our large-area surveys do not account for the Poisson nature of 
the detection. As we restrict the samples from the large-area sur¬ 
veys to highly significant detections, this simplification will have 
a minimal effect on our XLF measurements. However, differences 
in the assumed spectral shape can have a significant impact on the 
estimate of a flux and thus the assumed sensitivity. Uncertainties 
in the spectral shape and the resulting differences in sensitivity are 
accounted for by our Bayesian methodology described below. 


2.4 Identification and masking of stars 

Bright stars can contaminate our photometry, leading to issues 
with counterpart identification and photometric redshift estimates 
in these regions. We have therefore masked out areas close to bright 
stars from all of our Chandra fields in a consistent manner. We 
searched for stars brighter th an V = 15 in the HST Guide Star Cat¬ 
alog 2.3 jLasker et al.ll200d) . We masked all areas within a radius, 
r, given by 

r = (16 - U) X 6 " (1) 


Combining our five Chandra fields provides a sample of over 4000 
soft band detections and over 2800 hard band detections from 
a total area of ~ 9 deg^. However, to accurately constrain the 


where V is the V -band magnitude from the Guide Star Catalog. We 
set a maximum masking radius of 40”. We have removed any X- 
ray sources within this radius from our samples (which can include 
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the star itself or a nearby source). We also excluded these regions 
when calculating the X-ray sensitivity and area curves. 

We have also identified stars with fainter magnitudes that are 
detected at X-ray wavelengths and removed them from our sam¬ 
ples. In the fields with Rainbow coverage (CDFS, CDFN, EGS) 
stars were ident ified by a range of c olour and morphology criteria, 
as described in iBarro et alj ( 1201 lall . In the COSMOS and Bootes 
fields we applied a single colour criteri on based on the re gion of 
colour-colour space occupied by stars in lllbert et al.l l l2009l) . 

R - [3.6] < 3.0 X (R-I)- 1.2 (2) 

where [3.6] is the magnitude in the IRAC 3.6p.m imaging, R is 
the magnitude in the Subaru r+ or NDWFS R filter, and / is the 
magnitude in the Subaru i+ or NDWFS I filter. For all fields, we 
also required that the X-ray sources exhibit a low X-ray-to-optical 
flux ratio, log /x//opt < — 1, where the ratio is calculated as 

log = log /o.5-2keV + (3^ 

This cut ensures that we do not exclude bright QSOs from our sam¬ 
ple that may satisfy the other stellar criteria. When a spectroscopic 
classification is available (see Section iTSl this over-rides our pho¬ 
tometric classification. Of the 50 spectroscopically classified stars 
in our five Chandra fields, 40 (80 per cent) were also identified by 
our photometric procedure. 


2.5 Spectroscopic redshifts 


All of our Chandra surveys have been the subject of intense spec¬ 
troscopic campaigns.These campaigns include large-scale follow¬ 
up of the general galaxy population, in addition to those directly 
targeting X-ray sources. 

In the CDF- S we first s earche d for spectroscopic redshifts in 
the catalogue of IXue et al.l ( l201lh . The majority of the spectro¬ 
scopic redshifts for our CDFS-4Ms survey are taken from this cata¬ 
logue, as well as some of our ECDFS survey (in the area that over¬ 
laps with the 4Ms data). We also searched for spect roscopic red¬ 
shifts from the Arizona CDFS Environment Survey dCooper et al.l 
|2012|) . t he spectroscopic sub-sam ple of sources from the MUSY C 
sample dCardamone et alj2010d) . and PRIMUS dCoil et al.ll201lh . 

In the CDF-N we used spectroscopic r edshifts from DEEP 3 
([Coogeretalj201^ as well as the s ur veys o f|Trouine_et_^ l 
Barger et alj ( |20()^ : [Redd^_etay__ (2006); Wirth et ^ 

I.I d2004h andlstei 


2008h: 


Cowie et^ d2004 ) and lSteidel et alj ( 2003 ). 


2004 : 


In the ECS we compiled spectroscopic redshifts from a num¬ 
ber of surveys including DEEP2, DEEP3, the Canada-France 
Redshift Survey a nd MMT follow-up of X-ray sources. See 
iNandra et all d20l4 and references therein for full details. 

In COSMOS we initially searched for spectrosc opic redshifts 
of X- ray sources in the C-COSMOS catalogue of Icivano et al.l 
d2012h . We also search for additional spect roscopic redshifts from 
the bright zCOSMOS su rvey catalogue dUillv et alj 1 20091) and 
PRIMUS dCoiI etaI.ll20Ilh . 


In Bootes we adopt s pectroscopic redshifts f rom the AGN and 
Galaxy Evolution Survey dKochanek et alj[2012l) . For all fields we 
also searched f or spectroscopic r edshifts from the Sloan Digital Sky 
Survey (SDSS: lYork et al.ll200(]h . 

We matched the spectroscopic catalogues to the secure coun¬ 
terparts of our X-ray sources using a 2” search radius, corrected for 
any overall astrometric offset, and repeated the matching with a 1” 
radius. We only adopt those spectroscopic redshifts that are flagged 
as high-quality, reliable redshifts in the original catalogues. Table 


[T] gives the number of X-ray sources in our hard- and soft-band 
selected samples with reliable spectroscopic redshifts. 


2.6 Photometric redshifts 


The levels of spectroscopic completeness vary over our Chandra 
surveys from ~ 35 per cent (in our ECDFS area) to ~ 75 per cent 
(in the Bootes field). For the remaining X-ray sources we must re¬ 
sort to photometric redshift estimates, which are determined by fit¬ 
ting a set of template spectra to the observed spectral energy distri¬ 
butions (SEDs) of our sources. Such redshifts can be highly uncer¬ 
tain, particularly when considering faint X-ray sources. A key ad¬ 
vantage of our Bayesian analysis (see Section[3 is that we are able 
to account for the uncertainties in photometric redshifts by adopting 
probability distributions for the redshift, p[z), rather than a single 
redshift estimate. We thus require that our photo- 2 : approach recov¬ 
ers a p{z) distribution that accurately reflects the uncertainties in 
our redshift estimates. 

We calculate photome tric redshifts using the EaZY photo -2 
code dBrammer et al.ll200§) . We use EaZY in two-template mode, 
allowing for combinations of a galaxy and AGN template. For 
the galaxy templates we adopt the “pegasel3” template set pro¬ 
vided with EaZY. This template set consists of 259 synthetic 
ga laxy templates draw n from the range of parameters described 
bv ICrazian et alj l l2006l) . Additional “st ar-forming and dusty” tem¬ 
plates are included by applying the ICalzetti et alj ) l2000l) red¬ 
dening law for a range of different extinctions to a subset of 
the galaxy template s. We adopt seven AGN templates from the 
ISalvato et alj l l2009l) template seQ, namely the Sey 1.8, Sey 2, 


Mrk231, pLTQSOl, pLQSOH, pLQSO and the S0-10.QSO2-90 
hy brid template. We also include the Type-2 “Torus” template from 
the Pdlettaet^ d2007h libraryl^, which is not included in the fi 


Salvato et al. ( l2009h set. We note that some of these templates 


nal 

will include host galaxy contributions (particularly the Sey 1.8, 
Sey 2, Mrk231 and S0-10_QS02-90 templates). This is not a ma¬ 
jor issue as we want to determine the distribution of possible red¬ 
shifts, rather than perform an accurate host-AGN decomposition. 
The two-template mode in EaZY allows any possible combination 
of one of our AGN templates and one of our galaxy templates and 
thus allows for a large amount of flexibility in the fitted template 
SEDs. 

Another advantage of the EaZY photo -2 code is that it allows 
for the inclusion of a “template error function”. This feature ac¬ 
counts for additional uncertainty in the template SED as a func¬ 
tion of (rest-frame) wavelength and thus allows for the fact that our 
template set may not accurately represent the true diversity of SED 
shapes. This uncertainty in the true range of template SEDs is par¬ 
ticularly useful as our observed SEDs extend into the UV and mid- 
IR, where the templates are poorly calibrated, especially for AGNs. 
In addition, as the optical emission from an (unobscured) AGN can 
vary on timescales of months-to-years, our observed SEDs may not 
be well-matched by a single underlying template. The template er¬ 
ror function can also allow for any overall calibration uncertainties 
in the diverse sets of photometric observations used to construct 
our observed SEDs. We derive a template error function on a field- 
by-field basis—to ensure that it represents the calibration uncer- 
tainties in a given dat a set—using the basic procedure laid out in 
iBrammer et al.l ( 1200 8l) . Eirst, we attempt to fit the observed SEDs 


^ http://www.mpe.mpg.de/~mara/PHOTOZ_XCOSMOS/ 

® http://www.iasf-milano.inaf.it/~polletta/templates/swire_templates.html 
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Figure 2. Comparison of photometiic redshifts (^phot) spectroscopic redshifts (^spec) for our five Chandra fields. We plot the mean of the p{z) distribution 
as the best estimate of the photometric redshifts (black circles); error bars indicate the 95 per cent central confidence interval. The dashed line indicates a 1:1 
relation, whereas the dotted lines correspond to Hz/[1 + 2:zspec) = di0.15 (sources where the best estimate of the photo- 2 : lies outside this range are 
flagged as outliers). In the legend for each panel we give the number of X-ray sources with reliable spectroscopic redshifts (nzspec), the accuracy based on 
the normalized median absolute deviation (o-nmad)? the fraction of outliers (/outlier) the fraction of catastrophic failures (/catastrophic)^ see text for 
details. We highlight catastrophic failures with a red triangle. Over four of our fields (CDFS, CDFN, EGS and COSMOS) we obtain a consistent accuracy of 
CTNMAD ~ 0.05, with 15 per cent outliers and ~ 5 per cent catastrophic failures. In the Bootes field—our largest area, shallowest field—we have poorer 
accuracy and a higher outlier fraction, reflecting the more limited photometric imaging in this field. While our photo -2 have a poorer accuracy (ct^mad) 
a higher outlier rate than in some prior works, we choose to adopt our estimates as they have been calculated in a consistent manner across all five of our 
Chandra fields and have representative errors that we can fully track via the p{z) in our Bayesian methodology (see text for further discussion). 


with our templates, fixing the redshift at the spectroscopic value, 
where available. We then calculate 


A/.' = 


(4) 


where Fj indicates the observed flux in a filter for source j, and 
Tj is the flux from the best-fitting template. We calculate A/j as 
function of rest-frame wavelength for all sources and filters and 
take the median of every 400 individual measurements across the 
rest-frame wavelength range. We subtract the median photometric 
error, in quadrature, to estimate the contribution from “template 
error” to the uncertainty as a function of rest-frame wavelength. 
The value of the template error is typically around 10 per cent (in 
flux) but varies between ~ 4 per cent and ~ 20 per cent depending 
on the wavelength and the data in a given field. 

In Figure we compare our photo- 2 : estimates for X-ray 
sources to secure spectroscopic redshifts across our 5 Chandra 
fields. The panels for the EGS and CDFS fields combine the deep 


survey areas (AEGIS-XD, CDFS-4Ms), where the best photome¬ 
try is available, with the larger-area shallow surveys (AEGIS-XW, 
ECDFS). We include all X-ray sources with a high-quality spec¬ 
troscopic redshift in these plots and do not apply any cuts based 
on the estimated quality of the photo- 2 . Thus, we include sources 
with extremely broad p{z) distributions, which are often flagged 
as unreliable and excluded when assessing the success of photo -2 
techniques. 

In each panel of Figure we present a number of summary 
statistics: 

(i) ctnmad: the accuracy based on the normalized me¬ 
dian absolute deviation between the best photo -2 and 
the spectroscopic value, defined as ctnmad = 1.48 x 
median (|2phot - 2:spec|/(l + ^spec)). 

(ii) /outlier: the fraction of outliers, defined as the fraction of 
sources where |2phot — 2:speci/(l 4- •^^spec ) > 0.15; 

(iii) /catastrophic: the catastrophic outlier rate, calculated as the 
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Figure 3. Redshift distribution of X-ray sources in our five Chandra fields 
with high-quality spectroscopic redshifts (black solid line) and those where 
we adopt photometric redshifts. The blue dashed line indicates the distribu¬ 
tion of the “best” photo -2 estimates (mean of the p{z) distribution), whereas 
the orange dotted line shows the distribution obtained by combining the in¬ 
dividual p(z) distributions. 


fraction of sources where less than 5 per cent of the integrated 

p{z) lies within -0.15 < {z — Zapec)/Zapec < +0.15. 

Over four of our fields (CDFS, CDFN, EGS and COSMOS) 
we obtain a consistent accuracy of ctnmad ~ 0.06, with ~ 15 per 
cent outliers. Our approach ensures we assign an appropriate un¬ 
certainty, traced by the p{z), to the bulk of our sources and only 
~ 5 per cent of sources are thus flagged as catastrophic failures. 
In the Bootes field—our largest area, shallowest field—we have 
poorer accuracy and a much higher outlier fraction, reflecting the 
more limited photometric imaging in this field. However, the catas¬ 
trophic outlier fraction is only ~ 1 per cent, indicating that this 
additional uncertainty is represented by our p{z) distributions. We 
also note that in this field we have the highest spectroscopic com¬ 
pleteness (~ 75 per cent) and so only use photometric redshifts for 
a relatively small fraction of the sources in our eventual analysis of 
the XLF; when we do resort to a photometric redshift we account 
for the large uncertainty in the redshift. 

In Figure we plot the distribution of redshift estimates for 
sources in our Chandra fields where we have a spectroscopic red¬ 
shift and those where we adopt the photometric redshift informa¬ 
tion. Generally the sources with photometric redshifts lie at higher 
redshifts, a consequence of spectroscopic follow-up programmes 
being biased towards optically bright sources. We also show the in¬ 
tegrated contribution from the p{z) of all the photo-z sources. The 
integrated p[z) is skewed towards lower redshifts. This skew is due 
to the possibility that many potential high-redshift sources could 
actually lie at lower redshifts, which is reflected by their p{z) and 
must be accounted for in measurements of the XLF. 

A small fraction of our X-ray sources (< 2 per cent) lack a 
multiwavelength counterpart, precluding a photometric redshift es¬ 
timate. We retain these sources in our analysis, ensuring complete¬ 
ness of our sample, but adopt a p[z) that is constant in log(l + z) 
over our allowed redshift range (0 < z < 10). This reflects our 
lack of a priori knowledge of the redshift; the X-ray flux infor¬ 
mation is retained and thus a posteriori (after folding the constant 


p(z) through the final XLF) there may be a preferred redshift so¬ 
lution. The lack of a multiwavelength counterpart could imply that 
a high redshift solution should be given higher a priori preference 
for such sources, thus our approach is conservative. There is also 
a possibility that these X-ray sources lack counterparts as they are 
spurious detections, corresponding to positive fluctuations in the 
background count rate. Our analysis accounts for the Poisson na¬ 
ture of the X-ray detection and thus allows for this possibility. 

Other estimates of photometric redshif ts are available for 
X-ray sources in many of our fields (e.g. iBarger et alj 1200 3l ; 
ICardamone et alj|2010al : IHsu et al.l|2014h . Many of these studies 
take additional steps to improve the quality of the photo-z esti- 
mates. These s teps can include optimizing the template set (e.g. 
ILuo et alJ20i^ , attempting to corre ct the observed photometry for 
variability (e.g. ISalvato et alj|2009h . or appl ying priors based o n 
the source morphology and X-ray flux (e.g. ISalvato et al.ll201lh . 
These studies often achieve a higher accuracy and lower outlier 
rate than our own photo-z analysis. However, these additional steps 
can lead to underestimates of the true uncertainties in the photo-z. 
Conversely, we retain a large set of possible templates to ensure we 
produce p{z) distributions that account for the large uncertainties in 
the redshift and template degeneracies. While we have higher out¬ 
lier rates, our fraction of catastrophic failures remains < 5 per cent, 
indicating that ourp(z) distributions are representing the uncertain¬ 
ties. The nominal errors (i.e. the 68 per cent confidence intervals) 
on our photo-z are also comparable to the residuals between our 
best photo-z estimate and the available spectroscopic redshifts, in 
contrast to most previou s work where errors may be under estimated 
by a factor 2 — 6 (e.g. lLuo et al .l2010l:lHsu et al.l2014ll . Further¬ 
more, we require the full p(z) distribution, which most prior studies 
do not provide. We thus choose to use our own photometric red¬ 
shifts: our poorer accuracy and higher outlier rates are accounted 
for and compensated by our Bayesian analysis that incorporates 
the full p{z) information. 


3 BAYESIAN METHODOLOGY 

In this paper we expand on the Bayesian methodology developed 
by A10, incorporating the distribution of X-ray absorption proper¬ 
ties and accounting for the effects on the inferred shape and evo¬ 
lution of the XLF. Our method also accounts for the uncertainty in 
the measured X-ray flux (due to photon counting statistics), uncer¬ 
tainties in the redshift (for sources with photometric redshifts or no 
counterparts), uncertainties in the X-ray spectral shape, and the re¬ 
sulting uncertainty in the X-ray luminosity for an individual source. 
Our methodology is described below. 

3.1 Probability distribution function for a single source 

For a single source in either our hard- or soft-band sample we 
can derive the probability distribution function for z. Lx and A^h 
based on our observed data for that source alone, which is given 
by p{z, Lx, A^h I Di) where Di indicates the observed data from 
source i in our sample. This function is normalized such that 

dz j dlog Lx j dlogNn p{z,Lx,Nn \ Di) = 1. (5) 

We can re-write the probability distribution function as 

p{z, Lx, Nh I Di) = p{z I di) p{Lx, Wh | z,Ti, bi) (6) 

where di indicates the multiwavelength data used to estimate the 
redshift and Ti and bi correspond to the X-ray data for this source: 
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Rest-frame energy (keV) 

Figure 4. Our assumed X-ray spectral model for an AGN. We assume the 
intrinsic X-ray spectrum is a power law with photon index P = 1.9 di 
0.2 (purple line). The observed spectrum is absorbed by the intervening 
column density, fixed at Wjj = 5 x 10^^ cm~^ for this example (red 
dashed line). We also allow for a fraction (/scatt ~ 2 per cent) of the 
inti'insic power-law to be scattered, unabsorbed, into the line-of-sight (blue 
dot-dashed line). In addition, we allow for a component from Compton- 
reflection from cold, optically thick matter (such as a torus or accretion disc) 
that leads to the characteristic hump at high energies (green dotted line). 
The black line indicates the total observed spectrum, while the grey region 
indicates the 95 per cent confidence interval on the spectrum, allowing for 
the range of possible spectral parameters (see Section im Table^. 


the total observed X-ray counts in the given band and the estimated 
background respectively. 

Our knowledge of 2 is based on either a spectroscopic redshift, 
in which case we assume p{z \ di) is described by a J-function at 
the spectroscopic value, or a photometric redshift, when p{z \ di) 
is given by the p{z) from our photometric redshift fitting described 
in Section above. For the small fraction of sources where we 
were unable to identify a multiwavelength counterpart—and thus 
have no redshift information—we adopt 2ip{z) distribution with a 
constant density in log(l + z) over 0 < 2 < 10. 

The observed X-ray data can be described by a Poisson pro¬ 
cess. Thus, the likelihood of observing Ti counts from a source is 
given by 

£{Ti\ci,bi)= + (7) 

-t i ' 

where a is the X-ray count rate from source i in the observed 
energy band and we have assumed that the expected background 
count rate, bi, is well determined. The uncertainty in a is fully de¬ 
scribed by the Poisson likelihood given in Equation[7] However, to 
convert from a count rate, a, to an estimate of Lx and A^h (given 
z) we must assume a model for the X-ray spectral shape and fold 
this model through the appropriate instrumental response. 

Figure |4] shows an example of our X-ray spectral model. We 
assume the intrinsic X-ray continuum is described by a power law 
with photon index F and a high energy cut-off (we fix the fold¬ 
ing energy at 300 keV, although this has a negligible impact on 
the lower energies we observe). The observed spectrum is attenu¬ 
ated along the line-of-sight by the intervening column density, A^h. 
We include both photoelectric absorption (using the wabs model 
in XSpec) and Compton-scattering (via cabs), which suppresses 
the continuum further for Compton-thick column densities. A frac¬ 


tion of this power-law, /scatt, is allowed to emerge as an unab¬ 
sorbed component that is thought to be scattered into the line-of- 
sight by ionized gas in the vicinity of the AGN. We also include a 
contribution due to Compton reflection from cold, optically thick 
matter, which gives rise to a characteristic “hump” in the X-ray 
spectrum at 30 keV. Such a reflection component is expected 
due to reprocessing of the primary X-ray emission by a surround- 
ing, dusty torus or an accretio n disc. We adopt the pexrav model 
jMagdziarz & Zdziarskil [l995h . which is based on Monte-Carlo 
simulations of Compton reflection from a cold, optically thick slab 
of material. We assume the incident power-law has the same shape 
as the directly transmitted component. We fix the inclination angle 
to 30° as a representative value and allow the intensity to be set by 
the normalization, R, relative to that expected from a slab subtend¬ 
ing a solid angle of 27r (which is allowed to vary between 0 and 
2 , spanning the extreme cases of no reflection up to an effective 
Ftt solid angle coverage). We absorb the reflection component by 
the same column density seen by the primary emission. This is a 
good assumption when the reflection arises from the accretion disc 
and also provides reasonable agreement with the shape and inten¬ 
sity of the reflection co mponent based on sophistica ted models of 
toroidal obscurers (e.g [Srightman & Nandrair2011ah . All compo¬ 
nents are subjec ted to Galactic absorption, with column densities 
determined from lPickev & Lockmanl dl9^ via the Heasoft Nh 
tool. In XSpec terminology, our model is described by 

wabs * [ (1 — constant) * zwabs * cabs * zpowerlw * zhighect 

+ constant * zpowerlw 

+ zwabs pexrav] (8) 


where the constant corresponds to the scattered fraction, /scatt. We 
note that more physically motivated models could be adopted to 
describe the X-ray spectrum, self-consistently modeling the emis¬ 
sion, reflection_andabsorption_due_to_the_accretiondisc or torus 
(e.g. Ross & Fabia t]|2005l : lBrightman & Nandrj|201 lall . However, 
ISuchner et al. (2014) found that more simplistic models such as 
ours are generally sufficient to reproduce the observed spectral 
shape of individual, distant AGNs, especially considering our anal¬ 
ysis uses broad-band fluxes rather than performing a detailed X-ray 
spectral analysis. 

To describe our X-ray spectral model requires three additional 
parameters—the photon index, T, the scattered fraction, /scatt, and 
the relative normalization of the reflection component, R —which 
we introduce as “nuisance” parameters (collectively designated by 
4) in our Bayesian analysis. For a given set of spectral parameters 
(^), z. Lx and Nh, we can determine the expected count rate, a, 
and thus link the probability distribution function for Lx, A^h and 
^ to the Poisson likelihood given in Equation|3above. Thus, 


p(Lx, A^h,^ I z,Ti,bi) oc C{Ti \ Ci{z,Ly^,NH,i),hi)-K(f) 

(9) 

where Ci{z, Lx, Nh, $,) is the expected count rate for a source with 
redshift z, luminosity Lx, absorption column A^h, and additional 
spectral parameters ^ based on our X-ray spectral model, folded 
through the appropriate instrumental response for source i. 

A priori, the spectral parameters ^ for a given source are not 
well known; 7r(^) denotes the prior distribution that we adopt for 
our spectral parameters, which describes the range of possible val¬ 
ues and thus encapsulates the uncertainty in X-ray spectral shape. 
We assume the photon index, T, is drawn from a Gaussian distribu¬ 
tion with a mean of 1.9 and standard deviation of 0.2 (correspond¬ 
ing to the observed distribution of intrin sic photon indices i n X- 
ray spectral studies of nearby AGNs, e.g. lNandra et al.|[2007h . We 
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Table 2. Priors on the spectral parameters for an individual X-ray AGN, 4- 


Parameter 

Prior type 

Prior specification 

r 

Gaussian 

1.9 ± 0 . 2 “ 

/scatt 

lognormal 

-1.73 ±0.8*’ 

R 

constant 

0-2 


“Based on observed distribution fromlNan^aet^fcOOVl). 
*’Based on observed distribution from l^ntererai]j2009ll . 


assume the scattered fraction, /scatt, is drawn from a lognormal 
distribution with mean of 2 per cent and scatter of 0.8 dex based 
on the observed distribution o f partial covering fa ctors of sources 
in the SwiftiBKI sample from lWinter et al.l ( l2009l) . For the reflec¬ 
tion strength, R, we assume a uniform distribution in the range 
0 < R < 2. We thus allow for a large uncertainty in the strength 
of the reflection component, which is reasonable to encapsulate un¬ 
certainties in the geometry of the accretion disc and/or torus with 
our simplified modeling of the reflection. Table summarizes this 
prior information. The coloured lines in Figure |4] show each com¬ 
ponent evaluated at the prior mean of the spectral parameters and 
the black line corresponds to the total observed spectrum for these 
values. The grey region indicates the 95 per cent confidence interval 
on the total observed spectrum, adopting our prior distributions for 
the spectral parameters. We note that the large uncertainties in the 
scattered fraction, /scatt, lead to large uncertainties in the spectrum 
at softer energies for moderately and heavily absorbed sources and 
thus leads to large uncertainties in the intrinsic Lx- 

To obtain the probability distribution function for Lx and Wh 
only (for a given 2 : and our observed X-ray data), we must marginal¬ 
ize over the nuisance spectral parameters. Thus, 

p{Lx, Nh \ z,Ti,bi) — J p{Lx, \ z,Ti,bi). (10) 

In Figurewe show an example of p{Lx,Nii \ z,Ti,bi) for a 
source detected in the hard band (left) and a source detected in the 
soft band (right) with Chandra. In both cases we fix 2 = 1.0 and as¬ 
sume Ti = 30 total observed counts with an expected background 
of bi = 5.0 counts. We restrict the possible column densities to 
20 < log A^H(cm“^) < 26. The shading and contours indicate the 
range of possible values of Lx and Nu for each source. Detection 
in a single, broad energy band does not place constraints on the 
value of Nu, but we are able to place constraints on the range of 
possible values of Lx and how this depends on A^h . For the hard 
band example, we can constrain Lx to within ~ 0.25 dex, pro¬ 
vided the absorption column is Wh < 10^® cm~^. If the column 
density is higher, then the same observed counts must correspond to 
a higher value of Lx - For a source detected in the soft band, absorp¬ 
tion effects are apparent for lower column densities (Wh > 10^^ 
cm“^). At column densities A^h ^ 2 x 10^® cm“^ (at 2 = 1.0) the 
observed flux in the 0.5-2keV band is dominated by the scattered 
component only; the intrinsic luminosity is poorly constrained but 
must be a factor ~ 30 greater than if the source had a lower Nu. 

We note that the effect of absorption on the observed counts 
in the hard or soft energy bands will vary significantly with red- 
shift. For example, at 2 ~ 3, where the observed 0.5-2 keV energy 
band probes rest-frame energies 2 — 8keV, the inferred Lx is 
only affected for column densities > 10^® cm~^. These redshift- 
dependent effects are fully accounted for in our calculation of the 
expected counts, Ci{z, Lx, Nh,^), using our X-ray spectral model. 

For the analysis in this paper, we treat each source in our hard¬ 


er soft-band samples as independent detections. Thus, the only X- 
ray data we use for an individual source is the detected counts (and 
the expected background), which as described above does not al¬ 
low us to constrain the value of Ah for an individual source. In¬ 
stead, our approach (see Sections [3.2l below) involves determining 
the overall distribution of X-ray luminosities and absorption col¬ 
umn densities (given by the XLAF) that correctly describes both 
our hard- and soft-band samples. Using this method, certain values 
of Lx and Ah may be disfavored a posteriori (after performing 
our full analysis of the overall sample). For example, higher val¬ 
ues of Lx may be disfavored a posteriori based on the shape of 
the XLF (that generally shows higher Lx sources should be rarer). 
Conversely, lower values of Ah may be disfavored if our overall 
XLAF requires a large fraction of absorbed sources. 

Nevertheless, it is worth noting that a more sophisticated anal¬ 
ysis could place constraints on the value of Ah (or indeed any of 
the other spectral parameters) for an individual source before per¬ 
forming fits to the overall samples. We could combine the infor¬ 
mation on the counts in both the hard and soft bands for an indi¬ 
vidual source to place constraints on the underly ing X-ray spec¬ 
trum, based on t his hardness ratio information (e.g. lHasingeJ2008l ; 
IXue et al.l|20l^ . Alternatively, the full X-ray spectrum could be 
extracted for a source and fitted with our spectral model or any 
other well-motivated model (e.g. iTozzi et ^l2006l ; [Buchner et al.l 
|2014) . However, accurately incorporating this information into our 
Bayesian analysis—and fully propagating the uncertainties in the 
estimated spectral parameters. Lx and Ah— is beyond the scope 
of this paper. 


3.2 The likelihood function for the overall sample 


Here, we describe how we combine the individual probability dis¬ 
tribution functions for each source in our sample with a given 
model of the XLAF to construct the likelihood function for our 

overall sample. _ 

Following AlO (see also lLore^l2004h . we assume our indi¬ 
vidual sources are effectively Poisson points drawn from a distribu¬ 
tion described by the XLAF and the overall sample selection func¬ 
tion. The expected number of sources in our hard-band sample is 
given by 

■^fields 

t/(Lx,«,AH |©)7r(0 ^ A,-(Lx,2,Ah,0 
j=i 

where t/(Lx, z, Ah | 0) is the XLAF—the differential co-moving 
number density of AGNs per logarithmic interval in Lx and Ah— 
and is described by a model with a given set of parameters, ©. We 
note that the XLF, (j>{Lx,z \ ©) = (i.e. the differen¬ 

tial co-moving number density of AGNs per logarithmic interval in 
Lx), can be recovered from the full XLAF by integrating over Ah . 
Thus, 

cf>{Lx,z \ &) ^ j dlog Ah t/)(Lx,2, Ah I ©). (12) 

We investigate and compare various different parametrizations of 
the XLAF (or the XLF directly), which are described in Sections|4] 
and^below. 

The distribution of possible spectral parameters is given by 
7r(^), the prior on F, /scatt and R that we describe in Section [TT] 
above. ^ is the differential comoving volume per unit area, as a 


.Vhard(©) = yd log Lx yd log Ah / 
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Figure 5. Example constraints on the rest-frame 2-lOkeV luminosity. Lx, and the line-of-sight absorption column density, Aiji, for sources detected in the 
observed hard band (2-7keV, left) and soft band (0.5-2keV, right) in a typical deep Chandra observation. In both cases, we assume 30 total counts detected 
in the observed energy band with an expected background of 5 counts, typical of our sample in the AEGIS-XD survey region. We fix the redshift at z = 1.0. 
The grey shading indicates the probability distribution function for Lx and Aijj, allowing for the Poisson uncertainty from the observed counts and assuming 
the X-ray spectral model shown in Figure|4] marginalized over the uncertainty in the “nuisance” spectral parameters (F, /scatt, and R). The solid and dashed 
lines indicate the 68.3 per cent and 95.4 per cent (i.e. 1 and 2 a equivalent) confidence intervals on the joint parameter space. Detection in a single band does 
not constrain the value of Wjj but does place constraints on the range of possible values of Lx- At low values of Aijj, the X-ray luminosity is constrained to 
within ~ 0.25 dex. However, for higher A^h values the X-ray luminosity must be correspondingly higher to produce the same number of observed counts. 
Absoiption affects the X-ray spectrum (and thus the constraints on Lx) for lower column densities in the soft band (JVh > 10^^ cm~^) than the hard band 
(A^h > 10^® cm-2) at z = 1.0. 


function of z. Aj{Lx, z, A^h, $) is the area for field j that is sen¬ 
sitive to a source of luminosity Lx, redshift z, absorption column 
A^h, and spectral parameters effectively giving the probability 
of a source entering our sample. To calculate Aj{Lx, z, Nu,^) 
we use our X-ray spectral model described in Section 13.11 above 
to convert from a given Lx, z, Nu and ^ to an expected count rate, 
a [z, Lx , A^h , $) We then convert the count rate to an “effective 
flux” using our standard conversion factors and determine the area 
sensitive to this flux from our sensitivity curves (see Section ItTI 
and Figure[TJ. This calculation is performed for each of our Aifleids 
fields and the total area is summed. 

The likelihood function may now be constructed from the 
product of the probability distribution functions for the individual 
sources in the s ample and the probability that no other sources were 
detected (AlO, lLoreddl2004l) . For a Poisson process, the probabil¬ 
ity of no detected sources (in the hard-band sample) for a given 
model XLAF is Thus, for the hard-band sample the 

likelihood is 


L {Dhard 

^hard 

n 


C|'\ p—Miard(®) 


© 


J d\og Lx J dlogNu J 


(13) 


where L!(Dhard | ©) is the likelihood of obtaining the observed 
data from the entire hard-band sample (ZJhard) for a given set 
of parameters (0) for the XLAF, and p{z, Lx, Nu \ Di) is the 
probability distribution function for z, Lx and A^h for an individ¬ 
ual hard-band source with observed data (Di), marginalized over 
the nuisance spectral parameters (^). The product is taken over all 
rt-hard sources in our hard-band sample. 

The likelihood function for the soft-band sample can be con¬ 
structed in a completely analogous manner to that described for 
the hard-band sample above. As we treat the hard- and soft-band 
samples as independent samples, the overall likelihood function for 
both samples is simply given by the product of the two likelihood 
functions. Thus, 

£(Z?hard, Osoft I 0) = /:(I>hard | 0) ^(-Daoft | 0) (14) 

where we assume the same underlying XLAF (with parameters 0) 
can describe the observed data from both samples. This approach is 
taken in Section|5]where we combine the two samples to determine 
the form of the XLAF that adequately describes our two samples. 


3.3 Including the contribution of normal galaxies 


p{z,Lx,Nii \Di)i>{Lx,z,Nii | 6 >) 


® We assume an appropriate Galactic hydrogen column density, 

for each survey. For the ROSAT surveys and AMSS we assume the median 

Agof the sample. 


So far we have implicitly assumed that all of the X-ray point 
sources in our two samples are actually due to the emission from 
AGNs, rather than an alternative process. We have already excluded 
stars from our sample, either by masking known bright stars or ac¬ 
cording to the criteria given in Section [2A1 However, we have so 
far made no attempt to exclude normal galaxies from our sam¬ 
ple, where the X-ray emission may be predominantly due to the 
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Table 3. Priors on the galaxy luminosity function parameters, uj. 


Parameter 

Prior type 

Prior specification 

log A (Mpc“® dex“^) 

lognormal 

-3.96 ± 0.25“ 

a 

Gaussian 

0.76 ±0.10“ 

logLo (erg s-l) 

lognormal 

41.15 ±0.14“ 

/3 

Gaussian 

2.7 ±0.5*' 

Za 

constant 

0.5-2.5 


“Based on measurements in iGeorsakakis et"^ i2QQ6h . converted to rest- 
frame 2-lOkeV luminosities and our parametrization of the galaxy lumi¬ 
nosity function ('Eauat ionllSl. _ 

^Weak prior, based on iGeorgakakis etalj 120071) . 


combined emission from X-ray binaries within the galaxy—with 
the t otal luminosity tracing the star-forma tion rate of the galaxy 
(e.g. lRanalli et al.ll2003l : lMineo et al.ll2014h —rather than due to an 
AGN. The total luminosity produced by a galaxy is usually much 
lower than could be produced by an AGN, even in galaxies with 
extremely high star-formation rates. Thus, a common solution is to 
simply cut any source with an X-ray luminosity below Lx < 10“^^ 
erg s“^, and not consider sources below this l uminosity any fur¬ 
ther in measurements of the XLF of AGNs (e.g. [Barger et al.l2005l : 
ISilverman et al.l[2008ll . Such a cut is fairly conservative as a high 
fraction of the sources below this Lx cut will be dominated by 
AGN emission rather than star formation. Including such sources 
can also improve constraints on the faint end of the XLF. 

In our Bayesian analysis we do not assign a single Lx to any 
individual source; we have a probability distribution describing the 
range of possible Lx. Thus, we are unable to apply a strict luminos¬ 
ity cut to identify and remove normal galaxies from our sample. In¬ 
stead, we account for the potential contamination of our AGN sam¬ 
ple by including the X-ray luminosity function of normal galaxies 
(hereafter, the “galaxy luminosity function” or GLF) directly in our 
analysis. We assume the GLF is described by a Schechter function. 


number of galaxies in our hard-band X-ray sample is given by 


A/gal,hard(^) — log Lx J ^ dz dT' 


( 17 ) 


■^fields 


(?)gal(Lx, 2 I 7r(rgal) ^ ^ Aj (Lx , Fgal) 
7=1 


An analogous expression gives the predicted number of galaxies in 
the soft-band sample, A/’gai,soft(ut). 

We assume the spectrum for a galaxy’s X-ray emission can be 
described by a single power law (subjected only to Galactic absorp¬ 
tion), and thus is described by a single parameter, the photon index 
Tgai- We adopt a Gaussian distributio n with mean 1.9 and standard 
deviation 0.2 (e.g. lYoung et aljr2012h as the prior, 7r(rgai), on the 
range of possible photon indices. 

Our likelihood function (for the hard-band sample) must then 
be modified to include the expected number of normal galaxies and 
allow for the possibility that an individual detection is associated 
with a galaxy. Thus, EquationlJ^is modified to 

L(Dhard I X 

^hard 

i = l 

where 

= J dlogLxJ dlogNu J ^dz (19) 

(1 - rii)piz, Lx, Vh I Di) t/)(Lx, z, Nu \ ©) 


corresponds to the probability that source i is an AGN, for a given 
realization of the XLF and Wh function, and 

P]®"" = J dlogLx J ^dz (20) 

Vi PgalC^.Lx I Di) (f>^^l{Lx,Z I Lj) 


'('gal (Lx, z 


I bj) d log Lx = 





(15) 


where (()gai(Lx, « | at) is the GLF for a given set of parameters, at. 
We allow for pure luminosity evolution of the GLF by allowing the 
characteristic luminosity, L,, to vary with redshift as 


logL.( 2 ) 


logLo -I- /31og(l -I- 2 ) if 2 < 2c 
logLo -I- ,01og(l -I- 2c) if 2 > 2c 


(16) 


where Lq is the characteristic luminosity at 2 = 0 and 2c is a cut¬ 
off redshift. The GLF function is thus described by five parameters, 
at = [A, a, Lc,,P, 2c]. We adopt fairly strong (Gaussian or lognor¬ 
mal) priors on the principal parameters A, a, and Lo that describe 
the GLF at 2 = 0, based on a previous stud y using a completely 
independent dataset jGeorgakakis et al]|2006L see Tabled- As the 
redshift evolution is less well constrained, we app ly a weak, Gaus¬ 
sian prior on /3 (based on iGeorgakakis et alj2007f) and adopt a con¬ 
stant prior on 2 c, allowing any value in the range 0.5 < 2c < 2.5. 
Ultimately, we allow our own data to constrain these parameters 
through our exploration of the posterior parameter space (see Sec- 
tion l3.4l below). 

We must now modify the likelihood function for our overall 
sample to include the contribution of normal galaxies. The expected 


corresponds to the probability that source i is a galaxy, for a given 
realization of the GLF. Pgai( 2 , Lx j Di) describes the probability 
distribution function for Lx and 2 under the assumption that source 
i is a galaxy and is marginalized over the single nuisance spectral 
parameter, Fgai. Thus, 

Pgal(2,Lx|A) oc p{z \ di) X (21) 

j dFgal jC.(Ti I c'(2,Lx,rgal),&i^ ^^(rgal) 

where c(( 2 , Lx, Tgai) is the expected counts based on our single 
power-law X-ray spectral model (for a galaxy with redshift 2 , rest- 
frame 2-10 keV X-ray luminosity Lx, and photon index Fgai) and 
C{Ti I c'( 2 . Lx, Tgai), bi) is the Poisson likelihood (cf. Equation 
|7]l. We retain the same redshift probability distribution function, 
p{z I di), for the galaxy as was assumed in the case of an AGN. 

The parameter pi in Equations [T^andl^represents any prior 
knowledge as to whether source i is a galaxy or an AGN. This prior 
knowledge could be based on some other available data; for ex¬ 
ample an optical spectrum, morphological information, or the mul¬ 
tiwavelength SED. However, for this work we simply assume no 
a priori preference for a galaxy or AGN; thus we set t;; = 0.5 
for all sources. Instead, we allow our prior knowledge of the GLE 
itself (the relatively strong priors that we set on the parameters, 
u>) to drive our a posteriori inferences on the probability that a 
source is an AGN or a galaxy. As the GLF drops rapidly above 
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Lx ~ erg s“^, any X-ray source with a higher luminosity 

is very unlikely to be associated with a galaxy. Due to their unab¬ 
sorbed spectrum, galaxies will just as easily be seen in the soft band 
as the hard band; indeed, our approach naturally accounts for the 
rise in the number of sources detected in the soft-band at the ve ry 
faintest fluxes (e.g. lGeorgakakis et alj2008l ; lLehmer et al.ll20l3) . 

As in Section 13.21 above, the final likelihood function 
for the full data from both the hard- and soft-band samples, 
£(Dhard, DsoSt \ 0, oj), is given by the product of the hard-band 
and soft-band likelihood functions. 


3.4 Parameter estimation and model comparison 


Having determined the likelihood function for our overall sample, 
our knowledge of the XLAF and GLF can be fully described by the 
posterior probability distribution function. 


p(©,a; I D,M) 


C{D I @,Lj) I M) 

p{D I M) 


( 22 ) 


where D indicates all of our observed data and we have introduced 
M to indicate the hypothesis that a particular model parametriza- 
tion for the XLAF and GLF describes our observed data. The model 
itself is described by the parameters ©, oj; our prior knowledge of 
these parameters for a given model is described by 


7r(©,ciJ I M) = 7r(© I Mxlaf) 7r(ciJ | Mgai) (23) 


where Mxlaf and Mgai indicate our models for the XLAF and 
GLF respectively. 

The denominator in Equation|22l p{D \ M), is known as the 
Bayesian evidence, Z. This factor is calculated by integrating the 
likelihood function over the prior parameter space: 


2 


p{D I M) 



dw C(D I @,u;) \ M). 


(24) 


The posterior probability for a particular model (given the data) 
can be calculated from the evidence and any prior knowledge of 
the probability that the model is correct: 

p{M I D) = p{D I M) 7r(M) (25) 


The ratio of the posterior probabilities for different models, known 
as the Bayes factor, can be used for model comparison; if neither 
model is favored a priori, then the Bayes factor is simply given 
by the ratio of the evidences, Z, of each model. We adopt this ap¬ 
proach to compare between different models for the XLAF. We re¬ 
port the difference in logaritmic evidence between two models. 


A In .2 = In .El — In Z 2 (26) 

p(Mi I D) 
p(M2 I D) 

where Z\ is the Bayesian evidence for model 1 (Mi) and Z 2 is 
the Bayesian evidence for model 2 (M 2 ). A difference in the loga¬ 
rithmic evidence of A In 2 > 4.6 (corresponding to posterior odds 
of 100:1) indicates very strong evidence in favour of the model 1, 
whereas A l n.2 < —4.6 indicates very strong evidence in favor 
of model 2 ( |jeffrevslll96lh . If —4.6 < Aln.E < 4.6 then, while 
we may still favour the model with the higher evidence, we cannot 
decisively rule out one model in favour of the other. 

This Bayesian model comparison approach fully incorporates 
uncertainties in the underlying parameters and naturally applies 
Occam’s Razor, favouring a simpler model with fewer free param¬ 
eters over a more complex one, unless the latter is required by the 



data (e.g. iKass & Raftervlfl?^ . Our different models for the XLF 
or the full XLAF are described in Sections |4] and [3 below (along 
with a description of the appropriate priors on the different param¬ 
eters). Our model for the GLF, described in Section [33] above, is 
kept the same throughout the analysis of this paper. 

To perform the integration in Equation |24] and thus calculate 
the evidence for a given model of the XLAE, we adopt the MULTI- 
Nest algorithm 1 Feroz et aP 20091). which extends the “nested 
sampling” approach of iskillin 'd (l2004l) to allow for efficient explo¬ 
ration of large, multimodal parameter spaces. We use the MULTI- 
Nest code v3.6 with 400 “live points” and an efficiency factor of 
0.3. The algorithm also provides estimates of the posterior proba¬ 
bility distribution for the parameters, p(©, u; | D, M), which we 
use to obtain our “best estimates” of the parameters for a given 
model. We choose to report the posterior mean for a parameter as 
our best estimate. To characterize the uncertainty in a parameter, 
we report the 68.3 per cent (i.e. Icr equivalent) central confidence 
interval, marginalized over all other parameters. 


3.5 Binned estimates 


To aid in the visualization of our results and the comparison be¬ 
tween different models, it is useful to produce binned estimates of 
the XLF (including the GLF at low luminosities) in fine bins of Lx 
and 2 for both the hard- and so ft-band samples. To achieve this, we 
adopt the Nahs/Nmii method dMivaii et al.ll200ll) . as expanded on 
in A10. This method compares the observed number of sources in 
a given Ly^-z bin to that predicted based on a model fit. The binned 
estimate is then given by 


Noha 

Wmdl 


(j>{Lb, 26 I ©) + Zb I d>) 


(27) 


where Lb and Zb are the luminosity and redshift of the center of 
the bin, (j>b is the binned estimate, and © and d> are our best (a 
posteriori) estimates of the parameters for the XLAF and GLF. 

The predicted number of sources in a bin, Wmdi, for a given 
model XLAF and GLF, can be calculated from Equations 1 121 and 
M restricting the integration to the appropriate range in Lx and 
2 . To calculate the observed number of sources. Noha, we must ac¬ 
count for the distribution of possible values of Lx (and in many 
cases 2 ) for an individual source. A single source can make a par¬ 
tial contribution to the effective observed number in various Lx-z 
bins. To calculate Noha, we sum the partial contributions of the in¬ 
dividual sources to each Lx -2 bin: 


EL piogihi r 

Noha = J2 dlogLx / 

.;_n d log 1/1.-, J 2] 


-/logLio 


dlogLx / dz p{Lx,z\ Di,&,Lb). 


(28) 


where logLio, logLhi, 210 and 2 hi indicates the limits of the bin 
and p{Lx, z \ is the a posteriori probability distribution 

for Lx and 2 for source i, given the data for that source, Di, and 
our best estimate of the overall XLAF and GLF (described by our 
best estimates of the parameters, ©, Cj). Thus, 


p{Lx,z I oc 

Pga.l{z,Lx I Di)(j)g^a.\{Lx, z \ ^ 


(29) 


-b J dlogNhi(piz,Lx,Nhi I Di)tP{Lx,z,Nii \ ©)^ 


where we have marginalized over the probability distribution for 
Ah for source i. We assign an error to each binned estimate based 
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on the approximate Poisson error in the effective number of ob¬ 
served sources in each bin. We only plot points for bins with 

Noha 1- 

Our binned estimates of the XLF calculated in this manner 
will depend on the underlying XLAF model. For example, an 
XLAF that requires a high fraction of heavily absorbed sources 
would increase the binned estimate of the XLF based on the soft- 
band sample, to account for absorbed sources that would by miss¬ 
ing from our sample. Furthermore, the individual binned estimates 
are not independent as a single source can enter multiple bins. Nev¬ 
ertheless, the binned estimates serve to visualize and compare how 
well the data from our soft- and hard-band samples are described 
by a particular model. 


4 X-RAY LUMINOSITY FUNCTIONS FROM THE 
HARD-BAND AND SOFT-BAND SAMPLES 

In this section we investigate and compare different models to de¬ 
scribe the evolution of the XLF, neglecting the effects of absorp¬ 
tion. For this initial study, we fit the hard-band and soft-band sam¬ 
ples separately and determine the model parameters for the XLF 
from each sample. We adopt a log-constant distribution for Ah 
over a limited range (20 < log Ah < 21 ) but otherwise use our 
Bayesian methodology described in Section above (accounting 
for the distributions of spectral parameters and including the GLF). 
Given the lack of absorption corrections, the measurements of the 
XLF presented in this section can be seen as tracing the “observed” 
luminosity, corrected to 2-lOkeV rest-frame values. In Sections [4. II 
to l4.3l below we describe our various parametrizations of the XLF. 
Section l4!4l summarizes our findings. 


4.1 Luminosity-dependent density evolution (LDDE) 


Following a number of studies of the XLF of AGNs based 
on bo t h soft and h a rd X- r ay selected sample s (e.g. iMivaii et al.l 
120001: lEbrero et al.l l2009l: lUeda et alj I20l4 . we first investi¬ 
gated luminosity-dependent density evolution (LDDE) model 
parametrizations. These models start with the assumption that the 
XLF at z = 0 can be described by a smoothly connected double 
power-law. 


(?)(Lx, z = 0) 


d^Lx,z = 0) 


d log Lx 



where 71 is the faint-end slope, 72 is the bright-end slope, L« is the 
characteristic break luminosity, and K is the overall normalization. 
The XLF is then modified by an evolution term, e{z, Lx), which is 
a function of both redshift and luminosity: 


()>(Lx, z) = 0(Lx, a = 0) e( 2 , Lx). (31) 


We adopt two different models for the evolution ter m. The first, 
hereaf ter referred to as the LDDEl model, is taken from IUeda et al.l 
( l2003h . and assumes that e(z, Lx) is a power-law function of (1 -f 
2 ), with different indices above and below a cutoff redshift, Zd, 
which is itself a function of 2 . Thus, 


e(z,Lx) 


(l+zY^ [z ^ 2 ci(Lx)] 

(1 + 2.i(Lx))=^ > ^cl(^^x)] 

(32) 


where 


Zci(Lx) 


Zcl [Lx ^ Lai] 

[ix<Lal]. 


(33) 


For our second paraterization (hereafter LDDE2), w e adop t 
the further refinements of LDDEl proposed bv lUeda et al.l ( 120140 . 
The evolutionary term is modified to allow for a stronger density 
evolution at the highest redshifts, above a second cutoff redshift, 
2 c 2 . Thus, 


e(z, Lx) 


J (l + Zci(Lx))« 
(1 + zci(Lx)y^ 


1+2 

1 + 221 ( 1 x ) 
l + Ze2(Lx) 

l-FScl(Lx) 


e2 

^2 


[Z < 2cl(I/x)] 

[Zcl{Lx) < Z < Zc2(Lx)] 
[Z > Zc2iLx)] ■ 

(34) 


where the cutoff redshifts are both functions of Lx, given by 


and 


Zcl (Lx) = 


2c2(Lx) = 


' ^cl 


[Lx ^ Lai 

^ ^cl 

(&)•■ 

[Lx < Lai 

Zc2 


[Lx ^ La2] 

2c*2 ( 

'Lxr= 

La2 1 

[Lx < La2] 


(35) 


(36) 


Furthermore, in the LDDE2 model the ei parameter is allowed to 
depend on luminosity, with the form 


ei(Lx) = et-|-;3i(logLx - logLp) (37) 

as proposed in iHasinger et alJ (l200.5h . 

Both of these LDDE models have a large number of param¬ 
eters: 9 parameters are needed for the LDDEl model, whereas 
the LDDE2 model has 15 parameters. Given the large parame- 
ter space, earlier studies often fix some of these parameters (e.g. 
lUeda et alJlioIl fix 6 of the parameters in the LDDE2 model). 
In our Bayesian analysis, we instead apply priors for all parame¬ 
ters but otherwise allow them to vary, thus allowing the increased 
parameter space to be accounted for and appropriately penalized 
when calculating the Bayesian evidence. The physical meaning of 
the parameters is somewhat obscure, thus we apply constant (or 
log-constant) priors over reasonable ranges given the form of the 
parametrization and the range of our data. Table |4]gives our best a 
posteriori estimates of the parameters based on either the hard-band 
sample or the soft-band sample, along with our prior limits and 
the Bayesian evidence for LDDE2 relative to the simpler LDDEl 
parametrization. We note that the GLF is also included in each of 
the fits, introducing additional parameters, although these are partly 
constrained by the priors given in Table 

As can be seen from Table [4] we have very strong evidence 
in favor of the LDDE2 model over the LDDEl model for both the 
hard-band and soft-band samples, indicating that the introduction 
of the additional parameters in LDDE2 is justified by the data. 


4.2 Luminosity and density evolution (LADE) 


The LDDE model substantially warps the shape of the XLF with 
redshift, in particular introducing a flattening of the faint-end slope 
at 2 ~ 1 (as well as a further turn up at the lowest luminosities). 
Whether the observational evidence supports such cha nges in the 
shape of the XLF has been a matter of debate (e.g. AlO, Aird et al.l 


l2008l:lFiore et al.l2012l:IUeda et al.l201 j:lBuchner et al.l2015h . Fur 

thermore, a physical interpretation of this complex parametrization 
is difficult. Therefore, in AlO we proposed a model where the shape 
of the XLF (the smoothly-connected double power-law given in 
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Table 4. Prior limits and best a posteriori estimates of parameters for the LDDEl and LDDE2 models for separate fits to the hard-band and soft-band samples 
(neglecting absorption effects). 


Parameter 

Lower limit 

Upper limit 

Hard band 

Soft band 




LDDEl 

LDDE2 

LDDEl 

LDDE2 

\ogK (Mpc“^) 

- 7.0 

- 3.0 

-5.63 ±0.07 

-5.72 ±0.07 

-5.87 ±0.05 

-5.97 ±0.05 

logL* (erg s“l) 

43.0 

46.0 

44.10 ±0.05 

44.09 ±0.05 

44.17 ±0.04 

44.18 ±0.03 

71 

- 1.0 

1.5 

0.72 ± 0.02 

0.73 ±0.02 

0.67 ±0.02 

0.79 ±0.01 

72 

1.5 

4.0 

2.26 ±0.07 

2.22 ±0.06 

2.37 ±0.06 

2.55 ±0.05 

ei (orej) 

2.0 

6.0 

3.97 ±0.17 

4.34 ±0.18 

3.67 ±0.09 

4.35 ±0.35 

Pi 

- 2.0 

2.0 


-0.19 ±0.09 


0.65 ±0.06 

logLp (erg s-l) 

43.0 

46.0 


44.48 ±0.44 


44.45 ± 0.53 

e 2 

- 5.0 

0.0 

-2.08 ±0.17 

-0.30 ±0.13 

-2.92 ±0.14 

-0.96 ±0.12 

ea 

- 10.0 

-5.0 


-7.33 ±0.62 


-7.84 ±0.51 

^cl 

0.4 

2.5 

2.02 ±0.09 

1.85 ±0.08 

2.27 ±0.09 

1.80 ±0.08 

^c2 

2.5 

3.5 


3.16 ±0.10 


3.09 ±0.11 

log Lai (erg s“l) 

43.0 

46.0 

44.71 ± 0.09 

44.78 ±0.07 

0.92 ±0.11 

44.92 ±0.12 

logLa 2 (erg s-l) 

43.0 

46.0 


44.46 ±0.17 


44.27 ±0.30 

ai 

- 1.0 

1.0 

0.20 ±0.01 

0.23 ±0.01 

0.18 ±0.01 

0.16 ±0.01 

Oi2 

- 1.0 

1.0 


0.12 ±0.02 


0.06 ±0.02 

AlnZ 



0.0 

±36.1 

0.0 

±74.8 


Eauationl30labove) is kept the same at all redshifts, but undergoes 
a shift in luminosity with redshift as well as an overall decrease in 
density. We refer to this model as Luminosity And Density Ev olu- 
tion (LADE, see also e.g. lYencho et al.ll2009llRoss et alj2013L for 
similar parametrizations). The luminosity evolution is achieved by 
allowing L* to change with redshift, 


log L* [z) = log Lo - log 


t Zc 

1 + z 


1 -t- Zc 

1 + Z 


(38) 


where pi and p 2 allow for a different evolution of above and 
below a transition redshift, Zc- The additional density evolution was 
introduced by allowing the normalization, K, to evolve as 


logK(z) = logKo+d(l + z). (39) 


In AlO, we found that the Bayesian evidence for the LADE model 
was comparable to the evidence for the LDDEl modeQ when con¬ 
sidering a hard-band selected sample of X-ray sources at z < 1.2 
and colour pre-selected samples at higher redshifts. 

Here, we have repeated the fitting of the LADE model with 
our updated samples and improved analysis (allowing for a distri¬ 
bution of spectral parameters and the contribution from the GLF, 
but assuming all sources are unabsorbed). Our best a posteriori 
estimates of the parameters for the LADE model (along with the 
prior constraints) are reported in Table[5l along with the Bayesian 
evidence, relative to the LDDEl model fit. Eor both our samples, 
we find strong evidence in favour of the LADE model compared 
to the LDDEl model, indicating that it provides a better descrip¬ 
tion of the XLF despite having (slightly) fewer free parameters. 
The luminosity evolution of the best-fit LADE model in this work 
is somewhat different to the AlO findings, with a slightly weaker 
luminosity evolution at low redshifts (pi ~ 4, compared to pi « 6 
in AlO) that continues to higher redshifts (Zc ~ 2, compared to 
Zc ~ 0.8 in AlO) and has a much stronger, negative luminosity 
evolution at higher redshifts (p 2 ~ — 2 , compared to p 2 ~ —0.2 in 


^ The evidence for LDDEl was higher but not by enough to decisively 
favour LDDEl over LADE. 


AlO). These differences indicate that the LADE provides a better 
description than LDDEl for the high-redshift evolution, which is 
probed by our updated samples. Nevertheless, the evidence for the 
LDDE2 model—which introduces an additional density evolution 
at the highest redshfits—is significantly stronger than for the LADE 
model, despite the larger number of free parameters and additional 
complexity of LDDE2. 


4.3 Flexible double-power law (FDPL) 

The LADE model has a number of issues that may limits its ability 
to accurately described the observed behaviour of the XLF. Firstly, 
the forms of the luminosity evolution (in L,) and the density evolu¬ 
tion (in K) are restricted to the specified functional form. Thus, the 
model is unable to reproduce any more complex evolutionary be¬ 
haviour across certain redshift ranges, even if required by the data 
(e.g. at the highest redshifts). Secondly, the LADE model strictly 
requires that the shape of the XLF, described by the faint-end and 
bright-end slopes ( 71 , 72 ), remains exactly the same at all redshifts. 

We therefore propose a new set of models to parametrize the 
evolution of the XLF, which we refer to as the Flexible Double 
Power-Law (FDPL) models. These models assume that the XLF 
can be described by a double power-law form (Equationl30b at any 
redshift with four parameters: K, L«, 71 and 72 . We then allow 
any of these parameters to evolve with redshift, parametrizing the 
evolution by a polynomial function of log(l -|- z). By comparing 
the Bayesian evidence for models with different order polynomi¬ 
als, we can determine whether our data require a more complicated 
density or luminosity evolution, or whether there is any evidence 
for changes in the faint-end or bright-end slopes with redshift. 

With the simplest form of the polynomial, it is difficult to place 
meaningful priors on the free parameters: the various polynomial 
coefficients. Instead, we use Chebyshev polynomials over the in¬ 
terval 0 < log(l + z) < Cmax, where Cmax = log(l + Zmax) and 
we set Zmax = 7. For a given order, n, the polynomial form is de¬ 
fined by the value of the parameter (e.g. log K) at each of the n + 1 
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Table 5. Prior limits and best a posteriori estimates of parameters for the LADE model for separate fits to the hard-band and soft-band samples (neglecting 
absolution effects). 


Parameter 

Lower limit Upper limit 

Hard band 

LADE 

Soft band 

LADE 

log K (Mpc~^) 

- 7.0 

- 3.0 

-4.03 ±0.08 

-4.28 ±0.05 

log L* (erg s“l) 

43.0 

46.0 

44.84 ± 0.05 

44.93 ± 0.03 

71 

- 1.0 

1.5 

0.48 ± 0.03 

0.44 ± 0.02 

72 

1.5 

4.0 

2.27 ±0.07 

2.18 ±0.04 

Pl 

3.0 

10.0 

3.87 ±0.17 

3.39 ±0.08 

P2 

- 4.0 

3.0 

-2.12 ±0.39 

-3.58 ±0.26 

Zc 

0.4 

3.0 

2.00 ±0.13 

2.31 ±0.07 

d 

- 1.5 

0.5 

-0.19 ±0.02 

-0.22 ±0.01 

Aln.Z 



±9.6 

±38.4 


Chebyshev nodes, i.e. at redshifts, Zk, where 

1 1 / 2k 1 \ 

log(l + Zk) = -Cmax + -Cmax COS — -——TT . (40) 

A A \ A[n + i) / 

The value of e.g. log K{z) is then given by a sum of Chebyshev 
polynomials, 


log K{z) = '^CjTj{x) ( 41 ) 

j=o 


where Tj (x) are the Chebyshev polynomials of the first kind and x 
is a re-scaled version of = log(l + z) over the interval [—1,1], 




The coefficients, Cj , are given by 


( 42 ) 


n+1 


n + 


Y log A'(2; = Zk)Tk{x = Xk) 


( 43 ) 


where log iT(2 = 2^) is the value of log K at 2^. 

This scheme may appear complicated but it can easily be un¬ 
derstood as choosing values for a parameter (e.g. log K) at a num¬ 
ber of redshifts, Zk, and then finding a polynomial function to inter¬ 
polate between them for all 2. We can thus set a prior on the value 
of the parameter at each of the Chebyshev nodes {zk)- The choice of 
Chebyshev polynomials minimizes Runge’s phenomenon and thus 
prevents the adoption of highly oscillatory solutions for the redshift 
dependence. 

The FDPL models could take a large number of possible 
forms, with a polynomial of arbitrary order describing the redshift 
dependence of each of the four double power-law parameters {K, 
I/», 7i, 72). Fully exploring all possible combinations is computa¬ 
tionally prohibitive. To find an appropriate description of the XLF 
we proceed as follows. First, we try to find the best parametriza- 
tion for the overall density evolution and luminosity evolution by 
allowing log K{z) and log L* (2) to each be described by a poly¬ 
nomial of up to fourth order, but with no redshift-dependence for 
71 and 72. We evaluate the Bayesian evidence for all possible com¬ 
binations of polynomials of up to fourth order for log7C(2) and 
log L« (2) , corresponding to 25 different models. We find that the 
model where log K (2) is described by a second-order polynomial 
and logL«(2) is described by a third-order polynomial has the 
highest Bayesian evidence for the hard-band sample (although we 
note that alternative combinations cannot be ruled out). We then fix 


the polynomial orders for log 7^(2) and logL*( 2 j 3 evaluate 
the evidence for models where log 71 (2) (the faint-end slope) is 
described by a polynomial of up to fourth order. We choose to de¬ 
scribe log 71 (2) as a polynomial function, rather than 71 directly, 
so that extrapolation of our model to hig h redshifts does not r esult 
in negative, likely unphysical slopes (see lHonkins et al]| 2007 h . We 
find strong evidence that 71(2) changes with redshift, with a de¬ 
pendence that is described by a first-order polynomial (i.e. a linear 
relation) in log(l -I- 2); the higher-order polynomial dependences 
have lower Bayesian evidence but cannot be definitively ruled out. 
Finally, we fix the form for log 71 (2) to a first-order polynomial 
and test models where log 72 (2) (the bright-end slope) is allowed 
to vary with redshift. We find strong Bayesian evidence favouring 
the model with no redshift-dependence for 72. 

Our final “best-fit” model thus consists of log iT(2) being de¬ 
scribed by a second-order polynomial, log L* (2) being described 
by a third-order polynomial, log 71 (2) being described by a first- 
order polynomial and 72 being constant with redshift. Table|^gives 
the best a posteriori estimates of the parameters, corresponding to 
the values of the parameters at each of the Chebyshev nodes. We 
also provide a simplified polynomial expression for the redshift de¬ 
pendence of each of the double power-law parameters. In addition, 
we give the Bayesian evidence for this FDPL model, relative to the 
evidence for the LDDEl model (from Table| 4 ll. We find strong ev¬ 
idence in favour of our final FDPL model compared to any of the 
previously considered models, including LDDE 2 . 

We then repeat the entire process with the soft-band sample. 
We find that the soft-band sample requires the same model form 
as the hard-band sample, consisting of a second-order polynomial 
for log 7^(2), a third-order polynomial for log L,{z), a linear func¬ 
tion for log 71 (2) and a constant 72 (although the values of these 
parameters are different). The results, including the Bayesian ev¬ 
idence relative to the LDDEl model fit to the soft-band sample, 
are given in Table [ 7 ] Again, we find strong evidence favouring the 
FDPL model over LDDEl, LDDE 2 , or the LADE model. 


4.4 Summary of results 

In this section we have explored a number of different model 
parametrizations that can be used to describe the XLF, measured 


® The form of the redshift-dependence is thus fixed for log K{z) and 
log Lt{z) but the parameters that describe it are still allowed to vary in 
subsequent fits. 
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Table 6. Prior limits and best a posteriori estimates of parameters for the FDPL model for the hai‘d-band sample (neglecting absolution effects). 


Parameter 

Lower limit 

Upper limit 

A: = 0 

Parameter value at node k 

k = l k = 2 k = 3 

log A:(2 = 2fc) (MpC-3) 

-7.0 

-3.0 

-6.19 ±0.14 

-4.44 ±0.05 

-4.87 ±0.09 

logL«(2 = 2fc) (erg s“l) 

43.0 

46.0 

44.45 ±0.17 

44.61 ±0.05 

44.04 ±0.04 43.57 ±0.08 

II 

Fl 

0.01 

1.5 

0.27 ±0.04 

0.58 ±0.02 


72(2 = Zfc) 

1.5 

4.0 

2.31 ±0.07 




A In .3 +45.0 

log K(z) = -5.13 + 4.73C - T.IOC^ 
logL*(2) = 43.53 + 1.23C + 3.35C^ - 4.08C® 
log 71 (2) = -0.17-0.51C 

72(2) = 2.31 where C = log(l + 2) 


Table 7. Prior limits and best a posteriori estimates of parameters for the FDPL model for the soft-band sample (neglecting absorption effects). 


Parameter 

Lower limit 

Upper limit 

k = 0 

Parameter value at node k 

k = l k = 2 k = 3 

logL:(2 = 2fc) (Mpc“3) 

-7.0 

-3.0 

-6.47 ±0.09 

-4.73 ± 0.03 

-5.20 ±0.07 

logL«(2 = 2fc) (erg s-l) 

43.0 

46.0 

44.18 ±0.15 

44.61 ±0.04 

44.12 ±0.03 43.81 ±0.05 

FI 

II 

F) 

0.01 

1.5 

0.25 ±0.02 

0.55 ±0.02 


II 

Fi 

1.5 

4.0 

2.34 ±0.05 




A\nZ +93.1 

log K{z) = -5.47 + 4.88C - 7.20C2 

log L* (2) = 43.81 - 0.27C + 6.82C2 - 6.94C3 

log 71 (2) = -0.18 - 0.55C 

72(2) = 2.34 where C = log(l + 2) 


using either our hard-band or soft-band sample (neglecting the ef¬ 
fects of absorption). The analysis of both samples leads to simi¬ 
lar conclusions. We find that there i s stron g evidence for the up¬ 
dated LDDE2 model of lUeda et al.l ( 12014) compared to the sim¬ 
pler LDDEl parametrization or the LADE model of AlO (where 
the shape of the XLE remains the same at all redshifts). We also 
introduce a more flexible parametrization (FDPL) that models the 
XLE as a double power-law but allows for arbitrary evolution in the 
overall luminosity or density, as well as allowing for changes in the 
overall shape. The final FDPL model requires 10 free parameters 
to describe the XLE. Our FDPL provides a simpler parametriza¬ 
tion that is nevertheless able to reproduce the overall evolution of 
the XLE. The Bayesian evidence strongly favours our FDPL model 
over LDDEl, LDDE2 and LADE. 

In Figure 1^ we show how the double power-law parameters 
change with redshift based on our final FDPL model for the hard- 
band sample (black line). We find that the evolution of the XLE is 
driven by a combination of a) an overall density evolution (traced 
by the log K parameter) that peaks at 2 ~ 2, b) a relatively mild lu¬ 
minosity evolution that shifts the overall XLE towards higher lumi¬ 
nosities at higher redshifts (which continues out to high redshifts), 
and c) a mild flattening of the faint-end slope of the XLE with in¬ 
creasing redshift. For comparison, we also show the redshift de¬ 
pendence for the parameters based on our LADE model fit to the 
hard-band sample (dashed red line, along with uncertainties shown 
by the orange hatched region). The LADE model is unable to re¬ 


produce the form of density evolution required by our FDPL model 
and does not allow for any flattening of the faint-end slope. 

Figure [7] presents our best-fit FDPL model at a range of red¬ 
shifts based on the hard-band (blue line) and soft-band (red line) 
samples, along with binned estimates using the method described 
in Section 133] While we find that the form of the evolution of the 
XLF (described by the FDPL model) is very similar for both the 
hard-band and soft-band samples. Figure [7] shows that there are 
significant differences between the XLFs from the two samples. At 
high luminosities (Lx ^ L*) the space densities based on the hard- 
band and soft-band samples are broadly in agreement but at lower 
luminosities the best-fit model and binned estimates from the soft- 
band sample fall significantly below the hard-band estimates. This 
discrepancy is most likely due to our neglect of absorption effects 
and motivates further investigation in SectionObelow. 


5 JOINT CONSTRAINTS ON THE X-RAY LUMINOSITY 
AND ABSORPTION DISTRIBUTION FUNCTION 

In this section we attempt to find a parametrization for the full X- 
ray luminosity and absorption-distribution function (XLAF) that— 
after accounting for sensitivity effects and selection biases—can si¬ 
multaneously describe both our hard-band and soft-band samples. 
For each individual detection, we adopt a log-constant prior dis¬ 
tribution for Nh in the range 20 < log Ah < 26 (reflecting our 
lack of a priori knowledge of the absorption) and the corresponding 
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Figure 6. Dependence of the parameters that specify the double power-law form of the XLF on redshift based on best fits to the hard-band sample (neglecting 
absorption effects) with our FDPL model (black solid line) and LADE model (red dashed line). The solid grey and hatched orange regions indicate the 99 per 
cent confidence interval on the parameters for FDPL and LADE respectively based on the posterior distributions of the model parameters from our Bayesian 
analysis. The LADE model as proposed by AlO is restricted to a monotonic overall density evolution (traced by the log K parameter) and does not allow for 
any flattening of the faint-end slope (71). Thus LADE is unable to produce the evolutionary behaviour described by the more flexible FDPL model, which 
provides a better description of our data and is thus favoured according to the Bayesian evidence. 


constraints on Lx (e.g. Figure|5]l. Below, we describe our model for 
the XLAF (Section rS.lI l and present our results (Section r5.2l l. 


5.1 Model for the XLAF 

Most previous studies adopt a parametrization for the overall XLF 
of all AGNs (or, in many cases, all Compton-thin AGNs) and sepa¬ 
rately model the distribution of Ah at a given Lx and 2 : (often de¬ 
scribe d as the “Ah function” e.g. lUeda et al.l200^ : lLa Franca et al] 
l2005h . The overall shape of the Ah function and any luminos¬ 
ity or redshift dependence is defined by /abs, the fraction of ab¬ 
sorbed AGNfl Previous studies have shown that /abs is strongly 
depe ndent on luminosity and may als o increase at higher redshifts 
(e.g. Hasingejl2008l : lUeda et ^l2014l but see also Id welly & Pa^ 
|200( ). These dependencies may be modeled directly using specific 
parametrizations to describe /abs as a function of Lx and z; (e.g. 
power-law functions). Further assumptions, such as setting limits 
on /abs at certain luminosity or redshift thresholds, are often re¬ 
quired. Additional parameters are also required to fully describe 
the overall distribution of Ah for a given value of /abs- Such mod¬ 
els can have a very large number of parameters and thus many of 
the parameters may be fixed base d on measuremen ts in the local 
Universe or prior expectations (e.g. lUeda et al.l2()I^ . The resulting 
model parametrizations often have sharp breaks and discontinuities 
at certain luminosities or redshifts that may not be physical. 

In this work, we therefore choose to take a different ap¬ 
proach that builds on the flexibility of our FDPL model described 
above. We allow the unabsorbed (20 < log Ah < 22) and ab¬ 
sorbed (22 < log Ah < 24) AGN populations to be described 
by independent luminosity functions—<?!>unabs(Lx, z; | ©unabs) and 
((>abs(Lx, z; I 0abs)—with independent sets of parameters describ¬ 
ing their double power-law shape that can evolve differently with 
redshift. We assume a fixed fraction of unabsorbed AGNs, / 21 - 22 , 


® Defined here as the fraction of AGNs with column densities 22 < 
log Ah < 24 relative to all AGNs with log Ah < 24 


have column densities in the range 21 < log Ah < 22. Like¬ 
wise, we assume that a fixed fraction of absorbed AGNs, / 23 - 24 , 
have column densities in the range 23 < log Ah < 24. We as¬ 
sume a log-constant distribution of Ah within each of these 1 dex 
wide ranges (at a given Lx and z). The addition of / 21-22 and 
/23 -24 ensures we have sufficient freedom to describe the over¬ 
all distribution of Ah within the unabsorbed or absorbed popula¬ 
tions. We also allow for a population of Compton-thick AGNs with 
24 < log Ah < 26, which for this work are assumed to com¬ 
pletely track the evolution of absorbed AGNs. Thus, the XLF of 
the Compton-thick AGNs is assumed to be simply a factor, /3cthick, 
times the XLF of absorbed AGNs. Our overall model of the XLAF, 
^/(Lx, z, Ah I ©), is thus given by 


iP{L^,z,Nh I 0) 


(1 / 2 I — 22 )<(^unabs(Lx, Z | 0unabs) 

/ 2 I —22 *(*unabs(Lx, 2 | 0unabs) 

(1 — /23-24)?iabs(Lx, Z | ^abs) 
/23-24 *(>abs(Lx,Z | Pabs) 

<^Clh i ck 0„t„(Lx, Z I Paba) 


[20 


log Ah 

< 21) 

[21 


log Ah 

< 22) 

[22 


log Ah 

< 23) 

[23 


log Ah 

< 24) 

[24 


log Ah 

< 26) 




(44) 


where © = {0unaba, 0abs, /21-22,/ 23 - 24 , ;5cthick} and thus en¬ 
capsulates all parameters required to describe the XLAF. The factor 
i in Equation |44] for Compton-thick AGNs (24 ^ log Ah < 26) 
represents our assumption that the column densities are evenly dis¬ 
tributed over the 2 dex wide bin in log Ah. 


It is worth noting that in our analysis both the hard- and soft- 
band samples will contain AGN from the unabsorbed, absorbed, 
and Compton-thick AGN populations, although the soft-band sam¬ 
ple is likely to contain a lower fraction of absorbed sources than 
the hard-band sample (due to the detection biases against absorbed 
sources at soft energies). We do not know if an individual hard¬ 
er soft-band detection corresponds to an unabsorbed, absorbed, or 
Compton-thick source; we are only able to make statements about 
the population as a whole, based on our model for the XLAF. 
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Figure 7. XLFs calculated using the hard-band and soft-band X-ray samples, neglecting absorption effects and fitting each band individually. The solid blue 
and solid red lines show our best-fit FDPL models for the AGN XLF based on the hard-band and soft-band samples respectively. The green dashed line shows 
the GLF (soft-band fit only for claiity) and the dashed blue and dashed red lines show the total luminosity function, including both AGNs and galaxies, which 
should be compared to the binned estimates (dark blue triangles and orange circles for the hard-band and soft-band samples respectively). The sub-panels 
show the residuals between the data and model in terms of the A^obs/-^mdl ratio. While the form of the evolution based on the two samples is very similar, 
there are significant discrepancies in the space densities below the break in the XLF across all redshifts. To reconcile these discrepancies we must include the 
distribution of absorption column densities in our modeling and the effects on the inferred luminosities and space densities for both samples. 
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Marginalizing over log A^h, we can recover the total XLF, 

r26 

cj,{Lx,z\&)= / dlogN^ij{Ly,,z,N^\&) 

J 20 

— 0unabs(^X,2: | 0unabs) (45) 

4“ 0abs(-Ox5^ I ^abs) 

4“ /5cthick0abs (^X 5 ^ | ^abs) 

We assume that ijiunabs and 0abs are each individually de¬ 
scribed by an FDPL model paramaterization where the four XLF 
parameters— K, L,, 71 , 72 —can all vary with redshift according 
to polynomial functions of log(l 4- «)• Fully exploring all possible 
combinations of polynomials of different orders for the redshift- 
dependence of all eight of these shape parameters is computation¬ 
ally prohibitive. Instead, we initially assume that i^unabs and c^abs 
are both described by the form found in Section l43l when analyz¬ 
ing the hard- and soft-band samples separately. Thus, we adopt a 
second-order polynomial for log if ( 2 ), a third-order polynomial 
for log ( 2 ), a first-order polynomial for log 71 ( 2 ), and a constant 
72 as our baseline model. We then increase or decrease (if possible) 
the polynomial order for each of the eight shape parameters in turn 
and evaluate the Bayesian evidence relative to our baseline model. 


5.2 Results 

In Table[ 8 ]we give the change in the logarithmic Bayesian evidence, 
A In Z, found when increasing or decreasing the polynomial or¬ 
der for the redshift dependence of each of the FDPL parameters 
in turn, relative to the evidence for our baseline model. In the ma¬ 
jority of cases, the altered model has a lower Bayesian evidence, 
thus favouring our baseline model. For the normalization (K) of 
both the unabsorbed and absorbed AGN XLF and the L* of the ab¬ 
sorbed AGN XLF, we can decisively rule out the simpler models 
(with lower order polynomials used to describe the redshift depen¬ 
dence) as A In 2 < —4.6. For the L« of the absorbed AGN, a 
simpler (second-order polynomial) model is weakly favoured by 
our data (A In .Z > 0), but the evidence for this simplification is 
not strong. In all other cases, the baseline model is very weakly 
favoured (A In .E < 0). 

We also tested a model with no change in the faint-end slope 
for both the unabsorbed and absorbed AGN XLFs (but otherwise 
retaining the baseline model). While such a model is weakly fa¬ 
vored over the baseline (AlnZ = 0.71), we are unable to choose 
decisively between these models based purely on our data. We 
therefore retain the baseline model, which has greater flexibility 
to describe the faint end of the XLF. Nevertheless, as we are unable 
to fully explore all possible combinations of different polynomials, 
our final model must be seen as one possible, but adequate, way 
of describing our observational data. The evolution of the faint-end 
slope of the XLF is discussed in more detail in Section[6j2] below. 

In Tablel^we give our best estimates of the XLAF parameters: 
the shape parameters for the unabsorbed and absorbed XLFs and 
their redshift dependence in terms of Chebyshev polynomials and 
a simplified form, as well as f21-22, f23-21 and / 3 cthick- Table [Tol 
provides our best a posteriori estimates of the parameters of the 
GLF, which are constrained along with the AGN XLAF. 

Figure [ 8 ] shows our model for the total AGN XLF (found by 
marginalizing over A^h, see Equation 146b compared to binned es¬ 
timates based on our hard- and soft-band samples. The binned es¬ 
timates from both the hard- and soft-band samples are consistent 
with the total XLF (with an additional contribution from the GLF 
at low luminosities), showing that our overall model of the XLAF is 


able to describe both of our samples (cf. Figure^ where absorption 
effects were neglected). 

Figure 13 compares our model XLFs of the unabsorbed (blue) 
and absorbed (red) AGNs at a range of redshifts, as well as show¬ 
ing the total XLF of the full population (black). The XLF of ab¬ 
sorbed AGNs generally has a lower break luminosity (L*), a higher 
normalization, and a steeper faint-end slope than the XLF of un¬ 
absorbed AGNs. Hence, the absorbed AGN population is domi¬ 
nant at low luminosities, whereas unabsorbed AGN are increas¬ 
ingly important at higher luminosities. Both the absorbed and unab¬ 
sorbed AGN XLFs undergo a strong luminosity evolution, shifting 
to higher luminosities from 2 ~ 0 out to 2 ~ 3. Both XLFs also 
evolve in overall density, with the absorbed AGNs evolving more 
strongly and thus making up a larger fraction of the total AGN pop¬ 
ulation at 2 = 1 , even at high luminosities, but dropping away more 
rapidly at higher redshifts. The evolution of the XLF shape parame¬ 
ters with redshift is presented in Figure[T0] The differing evolution 
of the absorbed and unabsorbed AGN XLFs and thus their relative 
contributions to the total space density of AGNs results in a com¬ 
plex evolution in the shape of the total XLF of AGNs (solid black 

line in Figure 13- 

In Figure [TT] we plot the fraction of absorbed AGNs (rel¬ 
ative to the total Compton-thin population), /abs, as a function 
of Lx at various redshifts (i.e. the ratio of the absorbed AGN 
XLF to the sum of the absorbed and unabsorbed XLFs). The ab¬ 
sorbed fraction roughly plateaus at high and low luminosities, with 
a fairly rapid transition in between due to the difference in the 
break luminosities of the XLFs. The strong luminosity evolution 
of both XLFs causes this transition to shift to higher luminosities at 
higher redshifts. The stronger overall density evolution of the ab¬ 
sorbed AGNs also leads to an increase in /abs at high luminosities 
(Lx > erg s“^), where absorbed AGNs increasingly dom¬ 
inate, although our constraints on the bright end of the absorbed 
AGN XLF are fairly poor. Our model reproduces a luminosity 
and redshift dependence of the absorbed fraction and is generally 
consistent with direct estimates based on X-ray s pectral classifica¬ 
tions (e.e. lUeda et aflliol^ : iBuchner et al.ll2015ll . We do not find 
any evidence fo r the decrease in /abs with decreasing luminosity, 
which is seen by Buchneret 'ZIboish and som e other studies (e.g. 
iBurlon et alj201ll : |Briehtman & Nandral2011bh . although this may 
reflect the limited diagnostic power of our method to distinguish 
absorbed and unabsorbed AGNs at the lowest luminosities. 

Our model indicates that the positive luminosity evolution of 
the unabsorbed AGN population slows down by 2 ~ 3 and ap¬ 
pears to start a negative luminosity evolution to higher redshifts 
(see central panel of FieurefTOt. At 2 > 3, the negative luminosity 
evolution of the unabsorbed AGN XLF combined with the density 
evolution of the absorbed AGN XLF leads to a change in the pat¬ 
tern of fabs. At 2 = 3, /abs is approximately constant as a function 
of Lxlij, whereas at higher redshifts the absorbed AGNs appear to 
dominate at the highest luminosities. These results hint at a com¬ 
plex and differing evolution of absorbed and unabsorbed AGNs in 
the early Universe. However, our observational constraints on the 
total XLF at 2 ~ 4 — 7 are generally fairly poor. Furthermore, 
our diagnostic power to determine the absorbed fraction of AGN— 
comparing the hard- and soft-band X-ray samples—is substantially 
reduced at these redshifts. The hard-band sample is small at these 


Such a lack of lum inosity dependen ce in the absorbed fraction at 2 > 3 
has also been found bv ivito et alj 120141) . albeit classifying column densities 
of Ah > 10 ^® cm“^ as absorbed. 
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Figure 8. XLF of AGNs across the full range of absorption (20 < log A^h < 26) calculated by fitting both the hard- and soft-band X-ray samples to our full 
X-ray luminosity and absorption-distribution function (XLAF) model described in Section [sTI fsolid black line). We also show the GLF (green dashed line) 
which is fitted simultaneously with the AGN XLAF. Binned estimates are based on either the hard-band (dai'k blue triangles) or soft-band (orange circles) 
X-ray samples and should be compared to the total XLF, including both AGNs and galaxies (black dashed line). The sub-panels show the residuals between 
the data and model in terms of the Aobs/-^mdl ratio. Our final model of the full XLAF of AGNs allows us to reconcile measurements based on the hard- and 
soft-band X-ray samples as ultimately coming from the same AGN population (but subject to different selection biases and incompleteness as a function of 
Nn, z and Lx)- Hence, both the hai'd- and soft-band estimates are consistent with the total XLF (cf. Figure]?). 
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Figure 9. Our best model for the total XLF of AGNs (black solid line) and the constituent XLFs of the unabsorbed (20 < log A^h < 22: blue short-dashed 
lines) and absorbed (22 < log A^h < 24: red long-dashed lines) AGNs evaluated at a number of redshifts, based on our simultaneous fitting of both the 
hard- and soft-band X-ray samples to our XLAF model. The shaded regions indicate the 99 per cent confidence interval based on the posterior distribution 
of the model parameters. The hatched grey regions (and dashed black lines) indicate luminosity ranges where we lack either hard- or soft-band sources and 
our constraints are driven by extrapolation of our functional form. At 2 ;< 2, unabsorbed AGNs tend to dominate at high luminosties, whereas absorbed 
AGNs—due to the lower L* of their XLF—dominate at lower luminosities. Both populations undergo a strong luminosity evolution (generally shifting to 
higher luminosities as redshift increases) and an overall density evolution, changing the ratio of absorbed to unabsorbed AGNs at different luminosities and 
resulting in a complex evolution in the shape of the total XLF. Extrapolating our model to the highest redshifts {z ~ 5) indicates that absorbed AGNs may 
start to dominate at the highest luminosities; however, our observational constraints are weak in this regime and thus this behaviour may be an artefact of our 
particular model pai'ametrization. 



Redshift {z) 
0.0 0.5 1.0 2.0 



Figure 10. Shape parameters (normalization K, characteristic break luminosity L*, and faint-end slope 71) for the XLFs of unabsorbed (blue short-dashed 
line) and absorbed (red long-dashed line) as a function of redshift, calculated by fitting both the hard- and soft-band X-ray samples to our full XLAF model. 
The shaded regions indicate the 99 per cent confidence interval based on the posterior distribution of the model parameters and thus may under-represent the 
uncertainty in the observations at a particular redshift. 
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Table 8. Bayesian evidence (relative to our baseline XLAF model) for increasing or decreasing the polynomial order for the redshift dependence of each 
paramter in turn in the independent FDPL models (for the unabsorbed and absorbed populations) adopted for our modeling of the XLAF. 


AGN population Parameter 

Order 

(increased) 

A In .E Order 

(decreased) 

AlnZ 

Unabsorbed 

logA'(Mpc ^) 

3 

-0.86 

1 

-22.74 


logL* (erg s“l) 

4 

-3.12 

2 

-10.26 


71 

2 

-3.24 

0 

-0.25 


72 

1 

-2.45 



Absorbed 

log A” (Mpc“^) 

3 

-4.18 

1 

-27.23 


logL» (erg s“l) 

4 

-3.54 

2 

-12.85 


71 

2 

-4.23 

0 

-2.15 


72 

1 

-2.96 




Table 9. Prior limits and best a posteriori estimates of parameters for our full XLAF model, consisting of independent XLFs (described by the FDPL model) 
for the unabsorbed and absorbed AGN populations, and a scalar times the absorbed AGN XLF for the Compton-thick population. 


Parameter 

Lower limit 

Upper limit 

fc = 0 

Parameter value at node k 

k = l k = 2 

fc = 3 



Unabsorbed AGNs (20 < log Ajj 

< 22 ) 



logK{z = Zk) (Mpc-3) 

- 7.0 

-3.0 

-6.17L0.il 

-4.81 ±0.04 

-5.03 ±0.06 


logL «(2 = 2 fc) (erg s-l) 

43.0 

46.0 

43.26 ±0.14 

44.49 ± 0.03 

44.15 ±0.03 

43.79 ±0.05 

II 

N! 

0.01 

1.5 

0.04 ±0.02 

0.25 ±0.04 



II 

1.5 

4.0 

2.32 ±0.04 




/ 21-22 

0.1 

0.9 

0.43 ±0.04 




log K{z) = -5.21 ± 0.72C 

- 0.26C^ 






logL,( 2 ) = 43.56 ± 0.32C ± 0.32C2 _ 0.12C3 





log 71 (z) = -0.44-0.52C 







72 ( 2 ) = 2.32 



where = log(l + z) 





Absorbed AGNs (22 < log Nh < 24) 



logL :(2 = Zk) (Mpc-3) 

-7.0 

-3.0 

-6.81 ±0.09 

-4.44 ±0.05 

-4.30 ±0.09 


logL «(2 = Zk) (erg s”!) 

43.0 

46.0 

44.96 ±0.13 

44.57 ±0.06 

43.85 ±0.06 

43.17 ±0.08 

II 

0.01 

1.5 

0.16 ±0.04 

0.43 ± 0.04 



72(2 = 2 fc) 

1.5 

4.0 

2.33 ±0.18 




/23-24 

O.I 

0.9 

0.70 ±0.02 





log K{z) = -4.48 + 0.76C - 0.37C^ 
logL,{z) = 43.06 + 0.98C - + o.Olf^ 

log 7 i(^) = -0.28 - 0.28C 
72 ( 2 ) = 2.33 


where f = log(l + z) 


Compton-thick AGNs (24 < log A^h < 26) 
l3cthiok 0.2 2.0 0.34 ±0.08 


redshifts and both the hard and soft bands probe high rest-frame 
energies and are thus only weakly affected by absorption. In Fig¬ 
ure [TT] we indicate the range of luminosities where we detect both 
hard and soft band sources by the solid black line and grey con¬ 
fidence region. The dashed line and hatched grey region indicate 
where our constraints on /abs are primarily determined by extrap¬ 
olation of our functional form and should be treated with caution. 
Constraining the evolution of the bright-end of the XLF—in partic¬ 
ular testing whether the absorbed XLF evolves less rapidly than the 
unabsorbed XLF at 2 > 3, leading to the high /abs at Lx ^ 10'^® 
erg indicated by our extrapolation—requires large area, hard- 


band selected samples of X-ray sources, w hich will be provided by 
the eROSITA mission iMerloni et al.ll2012h . 


Nonetheless, our model is able to reproduce the overall evo¬ 
lution of the XLF and by accounting for absorption can recon¬ 
cile our hard- and soft-band X-ray samples. Furthermore, we natu¬ 
rally recover the luminosity dependence of the absorbed fraction at 
2 0 — 2 and its evolution. 
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Figure 11. Fraction of absorbed AGNs, /^bs (defined here as the fraction of AGNs with 22 < log < 24 relative to all Compton-thin, log Aljj < 24, 
AGNs) based on our best-fit model of the XLAF evaluated at a number of redshifts (solid black lines). The grey regions indicate the 99 per cent confidence 
interval based on the posterior distribution of the model parameters. The hatched grey regions (and dashed black lines) indicate luminosity ranges where we 
lack hard-band sources and thus our constraints on /^bs ate poor and driven by extrapolation of our functional form. The model evaluat ed at z = 0 is shown 
by the dotted line in all panel s. We also show direct estimates of the absorbed fraction as a function of luminosity from IUeda et al.l ^20141 green circles, for the 
indicated redshift ranges) and iBuchner et al.l 1201.4 orange shaded regions indicate their 90 per cent confidence intervals for bins that approximately span each 
redshift). Our flexible model reproduces a luminosity dependence of /abs. which plateaus at low and high luminosities, with a rapid transition in between. 
The position of this transition shifts to higher l uminosities at hig her r edshifts due to the lu minosity evolution of both the unabsorbed and absorbed AGN 
XLFs. This pattern is roughly consistent with the lUeda et al.l 120141) and iBuchner et alj 12014) measurements, although our estimates are slightly higher at low 
luminosities. At 2 = 3 we find that the absorbed fraction is roughly constant as a function of luminosity. At higher redshifts (2 = 5), absorbed AGNs may 
dominate at high luminosities, although this behaviour is based on extrapolation of our model XLFs to these high redshifts where we lack hard-band sources 
and are thus unable to directly constrain /abs- 


Table 10. Best a posteriori estimates of parameters describing the galaxy 
luminosity function. 


Parameter 

Value 

log A (Mpc“® dex“^) 

- 3.59 ±0.08 

a 

0.81 ± 0.03 

logLo (erg s-l) 

41.12 ± 0.07 

/3 

2.66 ± 0.24 

2 :c 

0.82 ± 0.10 


N.B. All parameters are partially constrained by 
our informative priors, specified in Table|3] 


6 DISCUSSSION 

6.1 Comparison with previous studies of the X-ray 
luminosity function 

Our work allows us to interpret the evolution of the total XLF of 
AGNs as due to the combination of the XLFs of unabsorbed and 
absorbed (as well as Compton-thick) AGNs. We find that both 
XLFs can be described by our FDPL model, albeit with differ¬ 


ent sets of parameters. Our approach differs from most prior stud¬ 
ies of the evolution of the XLF, which instead attempt to find a 
model that directly describes the total XLF, either using an LDDE 
para metrization or some form of luminosity and density evo lution 
(e.g. lEbrero et al.ll2009l : r^ncho et aklbOOgljAird et alj2010l) . 

In Eigure [T^ we compare the total XLF based on our fi¬ 
nal model (black solid line) to models from three recent studies: 
the LADE model from A10; the LDDE2 model put forward by 
lUeda et ^ 12014 herea fter U14); and the L DDE parametrization 
used in recent work bv iMivaii et alJ 120L5 . herea fter M15). We 
also compare estimates from iBuchner et af I 2 OI 5 I . herafter B15) 
who used a non-parametric method to estimate the space density of 
AGNs in fixed bins of Lx, and 2 with only simple smoothness 
assumptions used to link values between bins. 

At 2 < 3 there is good agreement between our work and the 
previous studies, although there are a number of important differ¬ 
ences. The bright end of our XLF is generally in good agreement 
with previous studies, although our updated work predicts a sub¬ 
stantially higher space density then the AlO model at both the low¬ 
est (2 = 0.1) and highest (2 > 2) redshifts. The B15 estimates 
are higher than our model at the bright end, especially at lower red- 
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shifts. However, this work lacked the very large-area (> 50 deg^) 
samples included in our study and other work. Thus, the B15 esti¬ 
mates at high luminosities and low redshifts may be driven by the 
smoothness assumptions inherent to their method. 

At 2 « 0.5 — 1 the faint end of our total XLF generally lies 
slightly above the AlO model. The differences between AlO and 
our updated study are likely due to a combination of absorption ef¬ 
fects (neglected in AlO), the limitations of the LADE model, and 
the smaller sample in the AlO work. The LDDE models from U14 
and Ml5 are in better agreement with our work, although there 
are slight differences in the shapes that may reflect differences in 
the parametrizations. For example, the M15 XLF is also somewhat 
higher than our model at the faintest luminosities at 2 ^ 0.5, which 
may be due to the LDDE parametrization requiring a steep faint- 
end slope (that is altered by the luminosity-dependent evolution). 
Our study probes to lower luminosities than Ml5, accounting for 
normal galaxies at low luminosities, and can thus rule out such a 
steep slope. We also predict a slightly higher space density than 
U14 or M15 at Lx ~ iq43-44 s“iand 2 = 0.5, and a slightly 
lower space density at moderate luminosities (Lx < 10'*'* erg s“*) 
at 2 = 3. The B15 estimates are generally in good agreement with 
our results at fainter luminosities out to 2 = 3, although our model 
tends to lie towards the lower end of their confidence intervals. 
There are small differences close to L* at 2 = 0.5 and 2 = 1.0 
whereby B15 predicts slightly lower space densities, which may 
also be driven by the assumed functional form and differences in 
the methodology. 

At Lx ~ lQi4,-45 gj.g g-l ^ _ g Q ig ^jg_ 

low the extrapolation of the AlO work or the B15 estimates but is 
consistent with the recent work of U14 and Ml5 (which did in¬ 
clude some sources at these redshifts and luminosities). We find a 
strong decline in the space density in of AGNs to high redshifts 
in this moderat e luminosity range, which is consist ent with previ¬ 
ous studies (e.g. ICivano et al.ll201ll : lvito et al.l2013ll . However, the 
high-redshift behaviour of our model may, in part, be driven by the 
extrapolation of the functional form that best describes the lower 
redshift XLE, where the bulk of our sample lies. We will present a 
focused study of the evolution of the XLE at high redshifts in future 
work (Georgakakis et al. in preparation). 

As an additional check, in EigurefT^we compare the extrap¬ 
olation of out model XLE to 2 = 0 with measurements based on 
local AGN populations. Our total XLE (left panel) does not undergo 
a sharp break, but flattens gradually below Lx ~ 10** erg due 
to the mixing of the unabsorbed and absorbed AGN XLEs. This 
mixing results in a relatively steep slope in t he Lx ]^q43-44 
erg luminosity range, consistent with the lUeda et alj l l201lh 
measurements shown by the black squares. The right panel com¬ 
pares direct measurements of the XLFs of una bsorbed (Wh < 10^^ 
cm“^) and absorbed (JVh > 10^^ cm~^) from iBurlon et alJ ^201 ih 
using the SwiftfBAT sample, where the populations were divided 
based on the results of X-ray spectral fitting for individual sources 
(cf. our statistical method). The extrapolation based on our base¬ 
line model for unabsorbed and absorbed AGNs is in good agree¬ 
ment with these measurements, which is reassuring given that our 
model is primarily constrained by much higher redshift sources 
and we do not directly measure A^h for individual sources. We 
note that the lowest two l uminosity bins for the unabsorbed AGNs 
from iBurlon et alj i201lh lie above our model (although only by 
1 — 2a). Thus. iBuidon et all ( 1201 ih found that the absorbed 
fraction decreases at both high and low lumin osities (see also 
iBrightman & Nandral201 Ibl i lBuchner et al.l20l3) . Our model does 


not generally predict such a behaviour in the absorbed fractiorF^ 
We find that /abs plateaus at ~ 85 per cent at low luminosities (see 
Figure fTTI l. In Figure [T^ (right panel) we also show our final best 
fit GLF. Normal galaxies have much lower lum inosities than the 
unabsorbed AGN found by [Burlon et al.l ( 1201 ih . Thus, the lack of 
a decrease in /abs at low luminosities does not appear to be due to 
us (statistically) misclassifying some fraction of our X-ray sources 
as galaxies rather than unabsorbed AGNs. 

In conclusion, our work is mostly consistent with recent esti¬ 
mates of the total XLF of AGNs at 2 « 0 — 4. However, our new 
model provides a simpler way of interpreting the evolution of the 
XLF and the absorbed fraction as due to the combination of the 
unabsorbed and absorbed AGN XLFs and their slightly different 
evolution. We discuss the possible physical origin of this behaviour 
in Section l63] below. 

6.2 Does the faint end of the XLF flatten at high redshifts? 

How the shape of the XLF evolves with redshift, and in particular 
whether the faint-end slope becomes progressively flatter at higher 
redshifts, remains a major question in studies o f the evolution o f 
the AGN population (e.g. iHopkins et alll2006l : lAird et al]|2010ll . 
In our study, we do not directly constrain the faint-end slope of 
the total XLF of AGNs - instead we describe the total XLF as a 
combination of the XLFs of the unabsorbed and absorbed AGNs 
and infer different faint-end slopes for each population (hereafter 
referred to as 7unabs and 7aba respectively). Our baseline model 
allows both 7unaba and 7ab8 to change with redshift and our results 
indicate that both slopes become flatter as redshift increases (see 
Figure □oj. At all redshifts, 7unaba is flatter than 7ab8. However, 
based purely on our data and analyses, we are unable to rule out a 
simpler model where neither 7unaba or 7ab8 change with redshift 
(hereafter the “no-change” model). We still require that 7unabs is 
flatter than y^ba in the no-change model. 

Despite these issues, we find that the faint-end slope of the 
total XLF does appear to change with redshift. This can be seen 
in Figure At 2 = 0.1 there is a gradual transition in the 
Lx ~ iq43-44 g-i range between the unabsorbed 

AGN XLF (dominant at high luminosities) and the absorbed AGN 
XLF (dominant at low luminosities) which results in in a relatively 
steep slope at these intermediate luminosities. It is only at fainter 
luminosities that the absorbed AGN XLF starts to dominate and 
a flatter overall slope is seerQ- At higher redshifts, the absorbed 
AGN XLF is increasingly dominant at moderate luminosities and 
thus a clearer break and a flatter faint-end slope is seen in the total 
XLF at Lx ~ iq42-44 ^-i py ~ 1, The underlying flatten¬ 

ing of 7ab8 in our baseline model contributes to this effect but the 
mixing of the populations appears to be more important. 

To further examine this effect, in FigurefTdlwe show estimates 
of the apparent shape of the total XLF as a function of redshift. We 
estimate double power-law shape parameters by taking our over¬ 
all model at various redshifts, creating fake data points in 0.25 dex 
wide redshift bins for 42.5 < log Lx < 46, and performing a 

We note that fitting the unabsorbed and absorbed XLFs as indpendent 
broken double power laws could reproduce such a pattern, as in iBurlon et alj 
1201 ih . thus the lack of such behaviour is not due to a limitation of our 
XLAF model. 

4^ Indeed, an advantage of our own study is that we probe down to very 
faint luminosities (Lx ~ 10^® erg s“4) at low redshifts using deep Chan¬ 
dra fields, carefully considering the contribution from normal galaxies, and 
improve the overall constraints on the faint end of the AGN XLF. 


© 2015 RAS, MNRAS 000. [TH3^ 






























26 


J. Aird et al. 



Figure 12. Our model for the total XLF of all Compton-thin (A^jj < 10^"* cm“^) AGNs (black solid line, grey region indicates 99 per cent confidence 
interval in model p arameters, dashed line and grey hatching indicate where we lack d ata and are extrapo lating our model) compared to prior model fits: the 
LADE model fromlAird et alj i2O10L blue long-dashed line); the LDDE2 model from|Ugda^yilJ i2014L green dot-dashed line); and the LDDE model from 
iMivaii et alj i2015l pink short-dashed line). We also show estimates from iBuchner et aIT]2oT^i ~who used a non-parametric method with simple smoothness 
assumptions to estimate the space density in fixed bins of Lx, N-n and z (orange hatched region indicate their 90 per cent confidence intervals for bins that 
approximately span the indicated redshift and include all Compton-thin AGNs). 


simple least squares fit to the standard double power law form. The 
solid black lines show the shape parameters of the total XLF for our 
baseline model, whereas the dashed green line shows the parame¬ 
ters for the no-change model. In both cases, the evolution of the 
total XLF is very similar - described by a combination of density 
evolution, positive luminosity evolution, and an apparent flattening 
of the faint-end slope as redshift increases. The apparent flattening 
is driven by the mixing of the constituent unabsorbed and absorbed 
AGN XLFs and their relative evolution. Both the baseline and no¬ 
change models predict a similar evolution of the total XLF, includ¬ 
ing this apparent flattening. We also show estimates of the shape 
parameters of the total XLF from Ml5, which they calculated by 
directly fitting the XLF with a double power-law in thin redshift 
shells. Their observed evolutionary pattern, including a strong flat¬ 
tening in the faint-end slope of the total XLF, is consistent with our 
findings. 

Additional evidence in support of our baseline model rather 
than the no-change scenario may be provided by the data shown in 
Figure [T3](right). These data are in good agreement with our base¬ 
line model, extrapolated to a = 0 (although the lowest luminosity 
bins for the absorbed XLF are systematically above our model). 
However, in the no-change model, 7abs and 7unabs are are driven 
towards the values required by our higher redshift data. Thus, the 
extrapolation of the no-c hange model to ^ = 0 does not provide a 
good agreement with the lBurlon et al.l (1201 ih measurements, lend¬ 
ing support to the baseline model. We n ote that in both ca ses our 
total XLF is in good agreement with the lUeda et alj 1 I2OIII) z = 0 


measurements (see Figure [T^ left). Incorporating these data sets 
into our Bayesian analysis remains beyond the scope of this paper. 

In conclusion, our own analysis requires that 7unabs is flatter 
than 7abs at any given redshift but does not provide definitive evi¬ 
dence that the individual slopes change with redshift. However, di¬ 
rect measurements of the XLFs of unabsorbed and absorbed AGNs 
at 2; = 0 indicate that redshift evolution of 7unabs and 7abs ft 
required when compared to our higher redshift study. Regardless 
of whether 7unabs and 7abs are changing, the apparent faint-end 
slope of the total XLF does flatten with increasing redshift. This 
behaviour is primarily driven by the changing mix of unabsorbed 
and absorbed AGNs. 


6.3 Evolution of the absorbed fraction of AGNs 

An increase in the fraction of absorbed AGNs with redshift (at 
a fixed luminosity) is now well establish ed (e.g. iLa Franca et al.l 
I2OO5I ; iHasingeilbOOgl ; iTreister et al.ll2009l) . Our model reproduces 
such a behaviour (see Figure [ 73 ] top). However, we are able to at¬ 
tribute the underlying cause of this evolution to the luminosity evo¬ 
lution the XLFs of both the unabsorbed and absorbed AGNs. This 
luminosity evolution shifts the XLFs of both populations to higher 
luminosities at higher redshifts, and thus the transition between the 
lower luminosity regime where absorbed AGNs dominate and the 
higher luminosity regime where unabsorbed AGNs dominate shifts 
to higher luminosities at higher redshifts (see Figure EB. Thus, 
close to the break luminosities of the XLFs (Lx ~ ]^q43.5-44.5 
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Figure 13. Left: Our best fit m odel for the total XLF of AGNs (black line) extrapolated to 2 = 0, compared to measurements of the local 2-10 keV XLF of 
AGNs from IUeda et alj ilOllh using MAXI (black squares). Right: Our best fit m odel for the XLF of unabsorbed (blue short-dashed line) and absorbed (red 
long-dashed line) AGNs extrapolated to 2 = 0, compared to measurements from iBurlon et alJ i2Ql ih from Swifi/BAT (converted to 2-10 keV luminosities 
assuming F = 1.9) where the populations were divided based on X-ray spectral fits (cyan diamonds and pink circles for unabsorbed and absorbed AGNs, 
respectively). The grey, blue and red regions indicate the 99 per cent confidence interval in our model extrapolation based on the posterior distribution of model 
parameters. Our model extrapolation, based on predominantly higher redshift samples and indirectly inferring the abso rption distribution, is gener ally in good 
agree ment with these local measurements. Our results indicate that the steep slope of the total XLF at 2 0 observed bv lueda ^ TiUToTil see also lBallantvnj 
|2Q14) may be due the gradual change i n the mix of absorbe d and unabsorbed AGNs. However, our model does not predict the increase in the unabsorbed AGN 
XLF at Lx $ 10^^ erg found bv lBurlQn~ 2 12 OI ih . which leads to a decreased absorbed fraction at low luminosities. 


Redshift (z) 
0.0 0.5 1.0 2.0 


5.0 


Redshift (z) 

i.O 0.5 1.0 2.0 5.0 0. 


Redshift (z) 

0 0.5 1.0 2.0 5.0 



Figure 14. The apparent double power-law shape parameters of the total XLF of AGN as a function of redshift. These estimates are based on a simple 
least-squares fit at various redshifts and in the luminosity range 42.5 < log Lx < 46.5 to the predictions of our model (which describes the XLF as the 
combination of independent XLFs of the unabsorbed and absorbed AGNs; see Figure FTol for the redshift dependence of the parameters of the individual XLFs 
that are directly constrained by our data). The solid black line indicates our baseline model where both 7unabs Tabs of fho underlying unabsorbed and 
absorbed AGN XLFs flatten with redshift, whereas the dashed green line uses our “no-change” model where Tunabs Tabs remain constant. The shaded 
regions represent 99 per cent confidence regions in the underlying models. In both cases, the apparent faint-end slope ( 71 ) undergoes a strong flattening 
with increasing redshift. The similarity in the parameters describing the total XLF for both models illustrates how distinguishing between these models is 
diffi cult without addition al information. Pink data points show estimates of the double power-law shape parameters of the total XLF in thin redshift shells 
from iMivaii et alj j2Q15h and reveal a similar evolutionary pattern. 
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Figure 15. Top: Absorbed fraction of Compton-thin AGNs (/abs) 
sus redshift based on our model, evaluated at three different luminosities. 
Shaded regions indicate 99 per cent confidence interval in the model param¬ 
eters. The absorbed fraction changes most rapidly at log Lx = 43.5 due to 
the evolution of L* for the unabsorbed and absorbed AGN XLFs. At high 
redshifts {z > 3) we have tentative evidence of a drop in /abs low and 
moderate luminosities and a rise at higher luminosities. Bottom: Fraction of 
Compton-thick AGNs (relative to all AGNs at a given luminosity) versus 
redshift from our model. For claiity, we omit the confidence regions for the 
low and high lumi nosities. We c ompa re to the estimat es of the Compton- 
thick f raction from iBurlonet green crosses), [Brightman & Ueifl 

b012L red circles'), and iBuchnereTal <20111 orange shaded region indi¬ 
cates 90 per cent confidence interval) that used X-ray spectral fitting to 
directly identify Compton-thick AGNs. While our estimate of /cThick ^t 
log Lx = 43.5 is systematically lower than prior works, we generally 
agree within the uncertainties. 

erg s“^) a strong evolution in the absorbed fraction is observed, 
whereas at higher and lower luminosities the evolution is weaker 
(see Figure[T5l top). 

Our results also suggest (albeit at low significance) that the 
absorbed fraction at high luminosities may increase slightly be¬ 
tween z ^ 0 and z 1 as the overall normalization of the ab¬ 
sorbed AGN XLF increases more rapidly than the normalization 
of the unabsorbed XLF. While such a redshift dependence is much 
weaker than the strong evolution around L*, it may explain why 
some studies have found substantial populations of obscured, lumi¬ 


nous AGNs (Tvpe-2 QSOs) at these high re dshifts (e.g. lSmsa et al.l 
l20ld : [^em et al.ll20ld : [Banerii et aDl2015h . The normalization of 
the absorbed AGN XLF peaks alz ~ 1 and declines towards higher 
redshifts. By 2 ~ 3 the continuing, strong luminosity evolution of 
the absorbed AGNs appears to take over, again leading to a high 
absorbed fraction at the highest luminosities (conversely the ab¬ 
sorbed fraction of lower luminosity AGNs may decline at z > 2). 
It should be noted, however, that our best-fitting model is primar¬ 
ily constrained by lower redshift data, where the bulk of our sam¬ 
ple lies. The high-redshift behaviour may simply reflect an extrap¬ 
olation of this model to higher redshifts where our observational 
constraints are poor. Nevertheless, our findings do hint at a poten¬ 
tially complex evolution in the absorbed fraction of high luminosity 
AGNs. 

6.4 The Nh distribution, the fraction of Compton-thick 
AGNs and their contribution to the cosmic X-ray 
background 

While we do not directly estimate absorption column densities for 
individual sources, our study does predict an overall distribution 
of Nu based on our statistical approach. In Figure [T^ we show 
our inferred distribution of Nu (the Nu function) evaluated at 
a = 0.05 and log Lx = 43.5. In the left panel, we compare to 
estimates of the Nh function from U14, based on di rect measure¬ 
ment s of Nh in the 5wi/I /BAT local AGN sample jTueller et al.l 
l2008tlidiikawa et al.l2012l) and corrected for absorption-dependent 
biases. Our inferred Nh function has a similar overall shape to the 
U14 measurements, dropping between logAln ~ 20 — 21 and 
log A^h ~ 22 — 23 and peaking again at log Nh « 23 — 24. 
Our absorbed fraction (/abs) is slightly lower at this specific lu¬ 
minosity and redshift, thus our Nh function predicts slightly more 
unabsorbed (log Nh = 20 — 22) sources and fewer absorbed 
(log Nh = 22—24) sources. Furthermore, our prediction of the A^h 
function in the Compton-thick regime (log Nh > 24) is lower than 
the U14 measurements and corresponds to an overall Compton- 
thick fraction (/cxhick) at z = 0.05 and log Lx = 43.5 of just 
11.4/g'g per cent. Even so, the differences between our model and 
the SwiftfBAT measurements from U14 are not significant. 

The right panel of Figure [T^ compares our results to the Nh 
functions assumed in various AGN population synthesis models, 
used to estimate the integrated cosmic X-ray background (CXB). 
While the overall shape of our Wh -function is similar, th e re are 
significant differe nces from prior stud ies. The iGilli et aP ( ItOOTL 
hereafter GOV) and iTreister et al.l ( l2009l . hereafter T09) models pre¬ 
dict fewer AGNs with Nh ~ ]^q 21-22 £^^-2^^ although neither 
our study nor the CXB models are particularly sensitive to the ex¬ 
act distribution in this range. More importantly, our A^H-function 
is significantly different in the Compton-thick regime. The normal¬ 
isation is significantly lower than the GOV or U14 models across 
the Nh = 10^°“^® cm~^ range, ruling out these distributions (and 
the correspondingly higher /cxhick) at the > 99 per cent confi¬ 
dence level. The T09 Wh -function, on the other hand, is higher 
than our model at Nh = 10^"*“^® cm~^ but assumes no sources 
with Nh > 10^® cm“^; thus, their overall /cxhick is consistent 
with our estimates where we assume a constant distribution up to 
Nh = 10^® cm“^. The exact shape of the A^h function in the 
Compton-thick regime, the relative number of Nh > 10^® cm“^ 
AGNs, and the precise X-ray spectral properties of such sources 
remain major uncertainties in estimates of /cxhick and could ac¬ 
count for the differences between our estimates and the prior stud¬ 
ies shown in Figure [16] (right). For example, if the most heavily 
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Figure 16. Our inferred distribution of absorption column densities (A^jj) evaluated at 2 = 0.05 and log Lx = 43.5 (black line, shaded r egion indicates 99 
per cent confidence interval in the model) and compared to direct measurements based on X-ray spectral fits of th e Swi^/BAT AGNs (left: lUeda et aljE014h 
and that assumed in previous CXB synthesis models (right: iGilli et aI]|2Q07l : lTreister et alj2009l ; IUeda et alj2014h . 


Compton-thick sources (log Nh > 25) are also deeply buried and 
do not exhibit the same level of soft, scattered em ission as less ab¬ 
sorbe d sources (as suggested by some studies e.g. lBrightman et Z1 
l2014h then our measurements would be biased to below the true 
value. 

Figure [Tsl (bottom) also compares our model for /cThi ck to 
estimates from lBurlon et alj i201lh , [Brightman & Uedal 1 2012h and 
B15 that used X-ray spectral fitting to directly identify Compton- 
thick AGNs in samples from Swifl/BAtT, Chandra, or XMM-Newton 
respectively. While our model is systematically lower than all of 
these measurements, it is generally consistent given the uncertain¬ 
ties in our model and the data. The B15 estimates are significantly 
higher at low (z < 0.2) and high {z > 3) redshifts, although these 
regimes are not well probed by their samples and thus the estimates 
may be partly driven by their priors. In the key 0.2 < 2 < 3 range, 
where the B15 estimates are most reliable, our estimates are statis¬ 
tically consistent, although the B15 confidence intervals do allow 
for the potential of a higher /cxhick. 

It is worth noting that our inferences regarding the Compton- 
thick population are indirect. Our estimates of the underlying dis¬ 
tribution of Nh rely on global comparisons of the hard- and soft- 
band samples and the changing sensitivities of these energy bands 
to different levels of absorption as a function of redshift. With 
this approach, our ability to distinguish between the effects of 
Compton-thick sources and heavily absorbed, yet Compton-thin 
(Nh ~ xq23-24 sources is limited, which could affect our 

estimates of /cxhick. We also assume that the XLF of Compton- 
thick AGNs traces the absorbed AGN XLF exactly (apart from a 
single scale factor) at all redshifts and assume a completely flat 


sues, our estimates of the Compton-thick fraction should be treated 
with caution. 

As an additional consistency check of our results, in FigurefTTl 
we use our model to estimate the overall contribution of AGNs to 
the CXB, retaining our assumed X-ray spectral model described in 
Section [TT] Our model prediction for the CXB is in good agree¬ 
ment with the measurements at energies > 10 keV that we do not 
directly probe with our source samples (we discuss the low energy 


CXB and the contribution of galaxies in Section 1631 below). De¬ 
spite our low Compton-thick fraction compared to U14 or GOV, 
we are still able to reproduce the observed peak in the CXB at 
~ 20 — 30 keV. At these energies, we find that Compton-thick 
AGNs make a small, but important contribution to the integrated 
CXB emission. 

Part of the reason for our successful reproduction of the CXB 
may be the comparatively strong reflection that we allow for in un¬ 
absorbed (log Nh < 22) and moderately absorbed (log Nh « 22— 
23) sources. The extent and strength of reflection in such sources, 
especially those with > L« luminosities and at 2 > 1 that make a 
large contribution to the CXB flux, is still a matter of debate (e.g. 
IVasudevan et alj2013l : lDel Moro et al.l2014l) . Indeed, the CXB can 
be reproduced with a wide variety of models that assume different 
X-ray spectral models for the constituent AGNs, as well as different 
XLFs, Nh distributi ons, and—perhaps most crucially—Compton- 
thick fractions (e.g . Draper & Ballantt^ I 2 OO 9 I : iBallantvne et al.l 
1 20 111 : lAkvlas et alll2012l) . While fully exploring CXB synthesis 
models for different spectral models and Compton-thick fractions 
is beyond the scope of this work. Figure [T^ does indicate that our 
relatively low Compton-thick fraction should not be ruled out. 


6.5 What drives the ahsorption-dependent evolution of the 
AGN population? 

The shape and evolution of the XLF of AGNs provides impor¬ 
tant constraints on the evolution of the underlying physical prop¬ 
erties of growing SMBFls, namely their masses and Eddington ra¬ 
tios (a tracer of the rate of accretion relative to the SMBH mass). 
Our study—revealing the independent XLFs of unabsorbed and ab¬ 
sorbed AGNs out to high redshifts—provides further insights into 
how the evolution of the AGN population relates to their absorption 
properties. 

We find that, at least to 2 ~ 2, the evolution of the unab¬ 
sorbed XLF and the evolution of the absorbed XLF are fairly sim¬ 
ilar. The bulk of the change in both XLFs is driven by moderately 
strong luminosity evolution, which could be attributed to a change 
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Figure 17. Predicted Cosmic X-ray Background (CXB) spectrum based 
on our final model of the XLAF and GLF (black solid line, grey region 
indicates 99 per cent confidence interval) compared to recent measure- 
ments from SwiftfXKl (orange circles at < 7 keV , iMoretti et alj |200^ 
and Swi/r/BAT (orange crosses at > 10 keV. lAiello et al.ll2008h . The blue 
doited, light-green dot-dashed, yellow triple-dot-dashed, and dark-red solid 
lines show the contributions from AGNs with different absorption column 
densities. The dashed purple line indicates the contribution of all Compton- 
thin (Alji < 10^^ cm~^) AGNs, while the solid purple line is the total 
AGN contribution including Compton-thick sources. The contribution of 
Compton-thick AGNs in our model to the peak of the CXB at ~ 30 keV is 
minimal but does bring the model closer to the observational data. The dark 
green long-dashed line indicates the contribution of normal galaxies to the 
CXB (based on our estimates of the GLF) and accounts for > 10 per cent 
of the total CXB flux in our model at < 3 keV. 

in the underlying distribution of SMBH masses, whereby the aver¬ 
age mass of an accreting SMBH is higher at earlier cosmic times. 
Alternatively, the luminosity evolution may reflect an increase in 
the average accretion rates with increasing redshift. Both XLFs also 
undergo an overall density evolution, whereby their normalizations 
increase between z = 0 and 2 ~ 1 then drop rapidly again at higher 
redshifts. This evolutionary pattern may reflect changes in the rate 
at which AGN activity is triggered within the galaxy population. Fi¬ 
nally, the faint-end slopes of both XLFs may also be changing with 
redshift. Such changes could indicate that the underlying distribu¬ 
tion of Eddington ratios is changing, with more rapidly accreting 
sources become relatively more predominant at higher redshifts, al¬ 
though recent studies find a fairl y constant, power-law distribution 
of Ed dington ratios to 2 ~ 1 — 2 jAird et al.l2012l : [Bongiorno et al.l 
|2012|) . The observed flattening of the total XLF may instead be 
more closely related to the changing mix of unabsorbed and ab¬ 
sorbed AGNs (see Section lb^ above). 

While the overall evolution of the XLFs of unabsorbed and 
absorbed AGNs follow similar patterns, the individual XLFs are 
significantly different at a given redshift. The most important dif¬ 
ferences (at least at 2 < 2) appear to be 1) the lower L* of the 
absorbed AGNs, and 2) the steeper faint-end slope of the absorbed 
AGN XLF. These effects combine with the overall evolution of the 
XLFs to produce the luminosity and redshift dependence of the ab¬ 
sorbed fraction seen in Figure fTTI 

The differences in L* may reflect the life cycle of AGN activ¬ 


ity. Luminous, unabsorbed AGNs could correspond to short periods 
in high accretion rate phases that have the power to eject material 
from the galactic nucleus. Exhaustion of the fuel supply then leads 
to a comparatively rapid fading of the AGN, that is traced by the 
faint end of the unabsorbed AGN XLF. Absorbed AGN, conversely, 
may correspond to longer periods when the SMBH grows at lower 
accretion rates but has little impact on the surrounding obscuring 
material. The two phases could correspond to separate AGN events 
with different triggering mechanisms (e.g. major mergers versus 
secular processes within a host galaxy). Changes in the extent and 
efficiency of different triggering mechanisms could then lead to the 
slight differences in the evolution of the unabsorbed and absorbed 
XLFs at 2 < 2. Alternatively, the unabsorbed and absorbed phases 
could reflect long-term variations related to a single AGN trigger¬ 
ing event. 

At 2 > 2, our results hint at a more complex evolutionary be¬ 
haviour. The luminosity evolution of the unabsorbed XLF appears 
to slow at 2 2 — 3 and then decline to higher redshifts, whereas 

the absorbed XLF continues to evolve positively in L* but the nor¬ 
malization (K) is rapidly declining (see Figure [TOt. This pattern 
suggests that absorbed growth phases continue to shift to higher 
accretion rates or SMBH masses at high redshifts, yet the trigger¬ 
ing of such events becomes increasingly rarer. Conversely, the most 
luminous, unabsorbed growth phases may be suppressed at the ear¬ 
liest times. However, our power to distinguish between the unab¬ 
sorbed and absorbed populations is limited at the these redshifts. 

The combined evolution of the unabsorbed and absorbed 
XLFs across our entire redshift range results in strong, luminosity- 
dependent evolution in the overall space densities of AGNs (see 
Figure [Tst. Higher luminosity AGNs peak, in terms of their space 
densities, at higher redshifts than lower luminosity AGNs. A simi¬ 
lar pattern has been seen in a number of previous works, most re¬ 
cently by M15 (data points in Figure fTSt . The minor differences be¬ 
tween the M15 estimates and our results are likely due to the differ¬ 
ent modelling of the underlying XLF and the correction for absorp- 
tionH This pattern in spa ce densities, described as “downsizi ng” 
of AGN luminosities (e.g. lUeda et al.ll2003l : lBarger et al.ll2005h . is 
often attributed purely to changes in the average SMBH mass with 
redshift. However, our constraints on the underlying evolution of 
the XLFs indicate that the observed pattern in the evolution of total 
space densities (Figure fT^ may be due to a combination of changes 
in SMBH mass, the distributions of Eddington ratios, and the life 
cycles of obscured and unobscured AGN phases. 

6.6 The contribution of normal galaxies to the observed 
luminosity function and the cosmic X-ray background 

An important advance in our work is the inclusion of normal. X-ray 
detected galaxies in our analysis of the XLE. Instead of applying a 
strict luminosity cut (e.g. Lx > 10"^^ erg s“^), we develop an ap¬ 
proach to statistically classify X-ray sources as galaxies or AGNs. 
We find that AGNs still dominate the space density of X-ray se¬ 
lected sources at Lx « 10"^^ erg at low redshifts. At 2 « 0.1 
the space density of galaxies is a factor ~ 1000 lower than the 
space density of AGNs at Lx = 10'*^ erg s“^. Galaxies only begin 
to make a substantial contribution (> 10 per cent in space density) 
at Lx < 2.5 X 10^^ erg s“^and start to dominate (> 50 per cent 
in space density) at Lx ;$ 2 x 10^° erg s“^. 

M15 use the Ajj-function from UI4 to correct for absorption effects (as 
a function of Lx and z), which differs from our results. 
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Figure 18. Total space density of AGNs for different ranges of X-ray lumi¬ 
nosity based on our model (coloured lines). Shaded regions indicate the 99 
per cent confidence interval in our model parameters. We see clear “down¬ 
sizing” in terms of AGN luminosity, whereby higher luminosity AGNs peak 
in terms of there space density at higher redshifts than low er luminosity 
source s. The data points are taken from the recent work of iMivaii et al J 
iZOlSli and reveal the same pattern; the small discrepancies with our model 
are most likely due to differences in the modelling of the XLF and the 
distribution. 


We find that the GLF undergoes strong luminosity evolution 
up to z ~ 0.8, whereby the characteristic luminosity increases as 
approximately (1 -|- zf’’^ (consistent with previ ous studies of the 
evolution of the GLF at X-ray wavelengths, e.g. Pteke^Q 2007|; 
iGeorgakakis et al] 1 20071 ; iTzmavaris & Georgantonoulosl l2008h . 
This evolution increases the contribution of galaxies at low lumi¬ 
nosities. At 2 ~ 1, the space density of galaxies corresponds to 
10 per cent of the AGN space density at Lx = 10^^ erg and 
galaxies dominate at Lx < 1.6 x 10“^^ erg s“^. Careful consider¬ 
ation of contamination by normal galaxies is thus vital in studies of 
faint AGN at these redshifts. 

At higher redshifts, our data are not deep enough to accurately 
probe the GLF and its evolution. Flowever, out model indicates that 
by z ~ 3 galaxies may dominate over AGNs at Lx = 10“*^ erg 
(see also Laird et ^120061 ; lAird et aDl2008h . Thus, galaxies may 
make a substantial contribution to the X-ray population at high red¬ 
shifts and should be carefully considered in future studies of the 
high-redshift XLF of AGNs. 

Galaxies also make a significant contribution to the CXB at 
low energies. Previous studies have shown that normal galaxies 
can account for ~ 4 — 20 per cent of the total CXB emission at 
1 — 6 keV energies, based on stacking of the X-ray data at 
the optical positions of known galaxies (e.g. IWorslev et aU 1 200a ; 
iHickox & Markevitd] |2007| ; IXue et al] l2012h . Flere, we can ac¬ 
count for the contribution of normal galaxies to the CXB based on 
our modeling of the resolved X-ray populations. Figure [771 shows 
that galaxies contribute > 10 per cent of the total CXB emission 
at energies < 3 keV. However, our estimate of the total CXB re¬ 
mains 10 — 20 per cent lower than the recent Swift/XKF mea¬ 
surements at < 3 keV, although measurements of the absolute flux 
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Figure 19. Bolometric luminosity density, a tracer of the total accretion 
density, for unabsorbed, absorbed (including Compton-thick) and all AGNs 
as a function redshift, based on our final model of the XLAF. Shaded regions 
indicate the 99 per cent confidence interval in our model parameters. The 
total rate of SMBH growth peaks at 2 : 2 and is dominated by absorbed 

and Compton-thick AGNs. 


of the CXB vary substantially between mis sions in this regime , 
probably due to cal ibration uncertainties (e.g. lBarcons et alfcood: 
livioretti et alJl2009l) . Indeed, our model of the total CXB is in ex¬ 
cellent agreemenTwitheMher measurements of the CXB flux with 
ASCA jGendreau et al.ll9^ over the full ~ 1 — 7 keV range. Thus, 
the significance of any discrepancy is uncertain. 


6.7 The accretion history of the Universe 

Our new measurements of the XLF of AGNs out to z 5 enable 
us to place improved constraints on the overall accretion history of 
the Universe. In Figure [79] we plot the total luminosity density of 
AGNs, 

y^dlogLboi<f>boi(Lboi,-z I 0) (46) 

where Lboi is the bolometric luminosity and (j>boi{Lhoi, z \ ©) is 
the bolometric luminosity function of AGNs, found by integrating 
our XLAF model over log A^h and converting the resulting XLF 
to bolometric values . We adopt the lumino sity-dependent bolomet¬ 
ric corrections from iHopkins et al] l l2007h . We find a clear peak in 
the luminosity de nsity at z = 1.76 ± 0 .05, consistent with prior 
studies (e.g. U14. lDelvecchio et al.ll2014l) . We also show the lumi¬ 
nosity density for the absorbed (including Compton-thick) and un¬ 
absorbed populations. At 2 > 0.1, the majority of SMBH growth is 
taking place in absorbed or Compton-thick AGNs. The proportion 
increases rapidly between z = 0 and the overall peak at a ~ 1 — 2. 
The contribution from unabsorbed AGNs increases more gradually 
and peaks at a slightly higher redshift (z = 1.89 ± 0.05). 

By integrating the AGN luminosity density across our red¬ 
shift range, we can estimat e the total SM BH mass density that has 
been built up by accretion iSoltanll 19821) . Adopting a radiative ef¬ 
ficiency of 77 = 0.1, we estimate a relic SMBH mass density of 
4.201 q j 4 ^ 10® Mq Mpc“® at z = 0, where the error repre¬ 
sents the uncertainty in our XLF model parameters. Of this mass 
density, ~ 30 per cent has been built up in unabsorbed growth 
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phases, ~ 50 per cent in absorbed AGNs, and a further ~ 20 per 
cent in Compton-thick sources. Our estimate is higher than in AlO 
(2.2±0.2 X 10® Mq Mpc“®), likely due to our improved account¬ 
ing for absorbed and Compton- thick sources, and is i n much better 
agreement with the estimate bv iMarconi et al.l l l2004h based on the 
velocity dispersions of local galaxies and the A4bh — cr relation 
(d.Olj ® X 10®; Al© Mpc“®). However, uncertainties in the bolo- 
metric correction and radiative efficiency dominate over our XLF 
model uncertainties. A higher radiative efficiency (e.g. if a large 
fraction of SMBHs have high spins) would reduce our estimate of 
the relic SMBH mass density. Alternatively, a large population of 
deeply buried, highly Compton-thick sources could bring our es¬ 
timate of the relic SMBH mass densit y more in line with r ecent 
revisions to the AIbh — cr relation (e.g. IComastri et alJl 2 oT^ . 

In Figurel20lwe compare estimates of th e total star formation 
rate ( SFR) density (from a recent review by iMadau & Dickinson! 
|2014) to our estimates of the SMBH accretion density, converted 
to Al© yr“^ Mpc~® using the same bolometric corrections and 
radiative efficiency assumptions described above and scaled up by 
an additional factor of 1500. We extrapolate our model to z > 5 
to enable comparisons with estimates of the SFR density in the 
first ~ 1 Gyr of cosmic time. The SFR density and SMBH ac¬ 
cretion rate both peak at 2 ~ 2. However, the total rate of SMBH 
growth rises slight more rapidly from 2 = 0 to a 2 , indicat¬ 
ing that the relative efficiency of SMBH growth relative to galaxy 
growth may change over this period. The SMBH accretion den¬ 
sity drops much more rapidly to higher redshifts than the SFR den¬ 
sity. Thus, galaxy growth may precede the build up of their central 
SMBHs in the early Universe, or the relative efficiency with which 
SMBHs form and grow may be substantially lower. However, our 
estimates are extrapolations at the highest redshifts. New and on¬ 
going surveys with Chandra (e.g. COSMOS-Legacy, P.I. Civano; 
CDFS-7Ms, RI. Brandt) may improve constraints in this regime. 
Nonetheless, obtaining direct and accurate measurements of the 
XLF of AGNs at the highest redshifts {z > 6 ) requires ~ 2 or¬ 
ders of magnitude improvement in survey speed and is one of the 
principal aims of ESA’s next large mission, the Athena X-ray ob¬ 
servatory jNandra et alj 20131 : 1 Aird et alj|2013bh . 


7 SUMMARY 

In this paper we present new measurements of the XLF of AGNs 
out to 2 ~ 5. In this section we first summarize our work and then 
present our conclusions. 

We compile large samples of uniformally selected X-ray 
sources in the hard (2-7 keV) and soft (0.5-2 keV) energy bands 
from a number of deep and wide Chandra surveys (CDFS, CDFN, 
EGS, COSMOS and Bootes), supplemented by wider area surveys 
from ASCA and ROSAT. We also compile multiwavelength pho¬ 
tometry across our fields that we use to robustly identify coun¬ 
terparts to our X-ray sources and calculate photometric redshifts. 
Next, we present a Bayesian methodology that allows us to incor¬ 
porate a range of X-ray spectral shapes and statistically account 
for the effects of absorption, as well as accounting for photometric 
redshift uncertainties, flux uncertainties, and the Eddington bias. 
We also introduce an approach to statistically account for the con¬ 
tamination of our samples by normal. X-ray detected galaxies. This 
allows us to place improved constraints on the XLE of AGNs down 
to much lower luminosities than prior studies. 

We present initial measurements of the XLF based on our 
hard-band and soft-band samples separately, neglecting absorption 


effects. We introduce a new, flexible model to describe the XLF of 
AGNs as an evolving double power-law, where the break luminos¬ 
ity (I/»), normalization (K) and faint-end slope ( 71 ) all evolve with 
redshift. There is significant evidence (based on Bayesian model 
comparison) favouring this model over luminosity-dependent den¬ 
sity evolution (LDDE) parametrizations of the XLF based on either 
the hard-band or soft-band samples. However, differences in the 
space densities at low luminosities based on the hard-band com¬ 
pared to the soft-band samples indicate that absorption has a sub¬ 
stantial effect and requires further consideration. 

To investigate the extent and distribution of absorption, we 
model the XLFs of unabsorbed (Ah < 10^^ cm“^) and absorbed 
(Ah > 10^^ cm“^) AGNs as independent, flexible double power- 
law functions. We use the combination of the hard- and soft-band 
samples to constrain our flexible model, seeking a model that can 
reconcile both samples assuming no a priori knowledge of the ex¬ 
tent of absorption for any individual source. We find that the ab¬ 
sorbed AGN population has an XLF that generally has a lower 
break luminosity, a higher normalization, and a steeper faint-end 
slope than the unabsorbed AGNs at all redshifts out to 2 ~ 2. Both 
the unabsorbed and absorbed AGN XLFs evolve in luminosity and 
density. Differences in the extent of the luminosity and density evo¬ 
lution lead to a comparatively complex evolution in the shape of the 
total XLF of AGNs. 

Our work provides new insight into the nature and evolution 
of the AGN population. Our main conclusions are as follows: 

1. Differences between the break luminosities (L*), normaliza¬ 
tions (A), and faint-end slopes ( 71 ) of the unabsorbed and ab¬ 
sorbed AGN XLFs explain the strong luminosity dependence 
in the fraction of absorbed AGNs at 2 < 2. Absorbed AGNs 
dominate at low luminosities, with the absorbed fraction falling 
rapidly as luminosity increases above a characteristic value. 

2. Both the unabsorbed and absorbed XLFs undergo strong lu¬ 
minosity evolution up to 2 ~ 2 , shifting the transition in the 
absorbed fraction to higher luminosities at higher redshifts. A 
strong evolution in the absorbed fraction is thus observed at 
Lx ~ 10 '* 3 . 5 - 44.5 gj.g g-i^ close to the break luminosities 
of the XLFs. At higher redshifts (2 > 2), the luminosity evo¬ 
lution of the unabsorbed AGN XLF slows and may become 
negative, although this behaviour may be driven by the model 
parametrization at lower redshifts. 

3. The normalization of the absorbed AGN XLF evolves more 
rapidly than the normalization of the unabsorbed AGN XLF, 
peaking at 2 ~ 1 and declining to higher redshifts. This leads 
to a mild increase in the absorbed fraction at high luminosities 
(Lx ~ 10 “^® erg s“^) between 2 ~ 0 and 2 ~ 1 . 

4. The total XLF of AGNs (combining unabsorbed and absorbed 
populations) undergoes a complicated evolution with redshift, 
including substantial flattening of the faint-end slope. This be¬ 
haviour is primarily driven by the changing mix of unabsorbed 
and absorbed AGNs. 

5. Our inferred distribution of Ah generally exhibits two peaks at 

Ah ~ iq 20-21 ^ although the exact 

shape depends on luminosity and redshift. We infer a relatively 
low fraction of Compton-thick AGNs (~ 10 — 20 per cent at 
Lx = lO'^®'® erg s“^). Nevertheless, our model reproduces 
the peak of the cosmic X-ray background at ~ 20 — 30 keV. 

6 . The majority of SMBH growth takes place in absorbed or 
Compton-thick AGNs at all redshifts 2 > 0.1, although the 
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Figure 20. Total star-formation rate (SFR) density compared to our estimate of the total SMBH a ccretion density. We show the best overall fit (dashed dark-red 
line) for the evolution of the SFR density from the recent review bv iMadau & DickinsonI i2014ll along with their compilation of measurements based on rest- 
frame ultraviolet (blue diamonds) and infrared (orange crosses) observations. Our estimate of the SMBH accretion density (solid black line) is scaled up by an 
arbitrary factor of 1500. Shaded regions indicate the 99 per cent confidence interval in our model parameters and may under-represent the overall unceifainties 
in the bolometric correction, radiative efficiency and observational data at a given redshift. The dashed black line and grey hatching indicate where we are 
substantially extrapolating our model outside the coverage of our data. Both galaxy and SMBH growth peak at 2 : ~ 2. However, the SMBH accretion density 
appears to evolve more rapidly, with a much stronger decline to higher redshifts. 


proportion increases rapidly between 2 = 0 and 2^2. The 
contribution from unabsorbed AGNs increases more gradually 
and peaks at a higher redshift (2 2). Unabsorbed phases of 

activity constitute a substantial fraction 30 per cent) of the 
overall accretion density. 

7. At low redshift (2 « 0.1) the space density of X-ray se¬ 
lected galaxies is negligible compared to AGNs at Lx = 10“*^ 
erg s“^. At 2 « 1, however, galaxies correspond to ~ 10 per 
cent of the AGN space density at Lx = 10“*^ erg and 
dominate at Lx < 1.6 x 10“^^ erg s“^. Normal galaxies con¬ 
tribute > 10 per cent of the total integrated CXB emission at 
energies < 3 keV. 

8 . The total AGN accretion density generally tracks the rise in 

the overall SFR density from 2 ~ 0 to 2 2, indicating 

that SMBFI growth and galaxy growth proceed together over 
this period. The accretion density appears to fall off substan¬ 
tially faster than the SFR density to higher redshifts, indicat¬ 
ing that galaxy growth may precede the build up of their cen¬ 
tral SMBFIs in the early Universe or the relative efficiency of 
SMBH growth may be lower. 

Our work provides new insight into the nature of the evolution 
of the AGN population and the underlying cause of the luminosity 
and redshift dependence of the absorbed fraction. However, to fully 
reveal the underlying physical processes that drive both unabsorbed 
and absorbed growth phases requires further study to connect the 
distributions of AGN accretion rates and their absorption properties 
to the key physical properties of the evolving galaxy population. 
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APPENDIX A: MULTI WAVELENGTH PHOTOMETRIC 
DATASETS 


In this appendix we provide details of the photometric datasets 
that are used to identify counterparts and calculate photometric 
redshifts for X-ray sources in our Chandra fields. The table be¬ 
low gives the name of all bandpasses and their effective wave¬ 
lengths. The prioritization in our likilihood ratio counterpart match¬ 
ing (where the band with priority 1 is used first and then the sub¬ 
sequent priority bands are used to identify additional, robust coun¬ 
terparts) is given in the third column. The final column gives the 
number of counterparts identified in each band and the correspond¬ 
ing fraction of the total X-ray sample. 

Band 

■^eff (^) priority A^cntrprt 


CDFS 

IRAC [3.6pm] 

35416.6 1 593 (95.6%) 

IRAC [4.5 pm] 

44826.2 

IRAC [5.8pm] 

56457.2 

IRAC [8.0pm] 

78264.8 

VIMOS U 

3711.4 

VIMOS R 

6414.0 2 7(1.1%) 

WFC3 F125W 

12425.8 

WFC3 F160W 

15324.7 3 0 (0.0%) 

ACS B 

4297.6 

ACS V 

5840.7 

ACS I 

7668.3 4 0 (0.0%) 

ACS Z 

9021.7 

ACS F814W 

7993.3 

ESO WFI U 

3414.4 

ESO WFI B 

4558.5 

ESO WFI V 

5355.9 

ESO WFI R 

6451.0 

ESO WFI I 

8568.1 5 5 (0.8%) 

CTIOz 

8966.4 

HAWK-I Y 

10194.5 

HAWK-I J 

12556.2 

HAWK-I K 

21423.2 

HAWK-I 1.061 pm 

10619.2 

SOFI J 

12444.9 

ISAAC J 

12459.9 

ISAAC H 

16453.2 

ISAAC K 

21594.9 6 0 (0.0%) 

SOFIA'S 

21591.4 

ISAAC Ks 

21594.9 

GALEX NUV 

2271.1 

GALEX FUV 

1528.1 

COMB017 B 

4554.4 

COMB017 V 

5358.1 

COMB017 R 

6432.4 

COMBOI7/ 

8523.6 

COMBO17F420M 

4177.7 

C0MB017F464M 

4616.8 
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Band 

Agfr (A) LR priority 

-A^cntrprt 

C0MB017F485M 

4858.5 



C0MB017F518M 

5188.1 



C0MB017F571M 

5716.1 



C0MB017 F604M 

6044.3 



C0MB017 F646M 

6450.9 



C0MB017F696M 

6958.9 



C0MB017F753M 

7530.5 



C0MB017F815M 

8157.2 



C0MB017F855M 

8556.4 



C0MB017F915M 

9140.0 



CDFN 

IRAC [3.6nm] 

35416.6 

1 

446 (90.1%) 

IRAC [4.5 urn] 

44826.2 



IRAC [5.8^m] 

56457.2 



IRAC [8.0^m] 

78264.8 



KPNOi7 

3566.7 



SUBARU B 

4366.4 



SUBARU V 

5429.6 



SUBARU R 

6486.4 



SUBARU I 

7952.0 

2 

13 (2.6%) 

SUBARU 2 ' 

9146.7 



ACS B 

4297.6 



ACS V 

5840.7 



ACS I 

7668.3 

3 

0 ( 0 . 0 %) 

ACS 2 

9021.7 



QUIRC HK 

18264.0 



CFHT J 

12518.7 



CAHA J 

12029.7 



CAHA iC 

22070.8 



CFHT JC 

21530.8 

4 

5 (1.0%) 

MOIRCS K 

21354.4 



CtAT.F.X FUV 

1528.1 



CtAT.F.X NUV 

2271.1 



WFC3 F105W 

10582.4 



WFC3 FI25W 

12425.8 



WFC3 FI60W 

15324.7 



EGS 

IRAC [3.6/.tm] 

35416.6 

1 

961 (92.5%) 

IRAC [4.5/im] 

44826.2 



IRAC [5.8/im] 

56457.2 



IRAC [ 8 |am] 

78264.8 



Subaru Re 

6486.4 

2 

51 (4.9%) 

ACS V 

5796.7 



ACS I 

8234.0 

3 

0 ( 0 . 0 %) 

Subaru Ks 

21354.4 

4 

1 ( 0 . 1 %) 

CFHTLS u* 

3805.6 



CFHTLS g' 

4833.7 



CFHTLS r' 

6234.1 



CFHTLS i' 

7659.1 

5 

1 ( 0 . 1 %) 

CFHTLS 2 ' 

8820.9 



MMT u' 

3604.1 



MMT g' 

4763.5 



MMT i' 

7770.5 

6 

1 ( 0 . 1 %) 

MMT 2 ' 

9030.9 



DEEPS 

4402.0 



DEEPS 

6595.1 

7 

10 ( 1 . 0 %) 

DEEP/ 

8118.7 



Palomar J 

12435.0 



Palomar K 5 

21353.0 

8 

4 (0.4%) 

CtALEX euv 

1528.1 




Band 

-^efl (A) LR priority 

-A^cntrprt 

GALEX NUV 

2271.1 


NICMOSEllOW 

10622.4 


NICMOS E160W 

15819.6 


CAHA J 

12029.7 


WEC3 E125W 

12425.8 


WEC3 E160W 

15324.7 


NEWFIRM Ji 

10441.4 


NEWFIRM J 2 

11930.0 


NEWFIRM J 3 

12764.0 


NEWFIRM Hi 

15585.4 


NEWFIRM H2 

17048.8 


NEWFIRM K 

21643.9 


COSMOS 

CFHT u* 

3805.0 


Subaru Bj 

4427.2 


Subaru Vj 

5454.9 


Subaru g'^ 

4728.6 


Subaru 

6249.1 


Subaru 

7646.0 1 

1338 ( 85.2%) 

Subaru 

9011.0 


IRAC [3.6gm] 

35416.6 2 

200 ( 12.7%) 

IRAC [4.5gm] 

44826.2 


IRAC [5.8gm] 

56457.2 


IRAC [ 8 gm] 

78264.8 


CFHT i* 

7582.6 


Subaru IB427 

4262.0 


Subaru IB464 

4633.7 


Subaru IB505 

5060.9 


Subaru IB574 

5762.8 


Subaru IB709 

7071.5 


Subaru IB827 

8242.5 


CFHT Ks 

21530.8 


UKIRT J 

12464.8 


Subaru IB484 

4847.6 


Subaru IB527 

5259.4 


Subaru IB624 

6230.9 


Subaru IB679 

6778.6 


Subaru IB738 

7359.5 


Subaru IB767 

7682.4 


GALEX NUV 

2271.1 


GALEX EUV 

1528.1 



Bootes 


NDWFS Bw 

4184.9 



NDWES S 

6501.8 



NDWES I 

8030.0 

1 

763 (91.7%) 

SDSSu 

3546.0 



SDSS g 

4669.6 



SDSSr 

6156.2 



SDSSi 

7471.6 

2 

6 (0.7%) 

SDSSz 

8917.4 



NDWFS K 

22078.4 

3 

9(1.1%) 

FLAMINGOS J 

12442.5 



FLAMINGOS Ks 

21483.9 



IRAC [3.6Atm] 

35416.6 

4 

40 (4.8%) 

IRAC [d.Sgm] 

44826.2 



IRAC [S.Sgm] 

56457.2 



IRAC [ 8 nm] 

78264.8 



GALEX FUV 

1528.1 



GALEX NUV 

2271.1 

5 

3 (0.4%) 
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